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ABSTRACT 


This  report  describes  improvements  which  have  been 
made  in  the  Airborne  Integrated  Reconnaissance  System  (AIRS) 
Performance  Model.  Specifically,  the  report  includes: 


. An  algorithm  for  analyzing  the  observability 
of  target  complexes. 

. Methodology  for  evaluating  the  effectiveness 
of  slewing  optical  sensors. 

. A technique  for' examining  why  sensors  perform 
as  they  do,  on  a micro-level. 

. Improvements  in  the  performance  models  of  in- 
dividual sensors  by  inclusion  of  additional 
equipment,  environmental,  and  operational 
parameters . 

. Improvements  in  the  method  of  data  presentation 
to  the  AIRS  user. 


In  addition,  the  prospective  AIRS  user  is  provided 
with  instructions  on  the  use  of  the  several  computer  programs 
which  make  up  the  overall  model. 
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SECTION  I 
INTRODUCTION 


1 . 1 BACKGROUND 


Under  NAVAIRDEVCEN  contract  number  N62269-68-C-0441  , 
Analytics  developed  a simulation  model  for  studying  the 
performance  of  future  tactical  airborne  reconnaissance  systems. 
This  simulation  model,  called  the  AIRS  Model  (for  Airborne 
Integrated  Reconnaissance  System) , was  documented  in  the 
following  report  to  NAVAIRDEVCEN: 

Analytics  Incorporated,  Final  Report  1004-1, 

Simulation  of  Advanced  Integrated  Reconnaissance 
Systems  (AIRS) , Volumes  I,  II,  and  III.  Submitted 
to  Naval  Air  Development  Center,  Johnsville, 
under  Contract  N62269-68-C-0441  , 4 April  1969. 

Under  a subsequent  contract,  Analytics  was  tasked 
to  make  improvements  in  the  AIRS  Model.  This  report  documents 
these  improvements. 

It  is  assumed  that  the  reader  has  a detailed  knowledge 
of  the  three  AIRS  report  volumes , this  report  cannot  be 
effectively  used  unless  this  is  the  case.  The  AIRS  report 
will  be  cited  frequently  throughout  this  report,  and  will  be 
referred  to  simply  as  "AIRS",  with,  of  course,  the  appropriate 
volume  and  paragraph  number. 
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It  would  be  good  to  redefine  here  certain  key 
concepts  of  the  AIRS  Model  (see  AIRS,  Volume  I,  paragraph  1.1). 
The  four  primary  measures  of  performance  for  tactical  air- 
borne reconnaissance  systems  used  by  AIRS  are: 

1)  Target  detectability. 

2)  Target  identif iability . 

3)  Localizability  of  the  target  (expressed  as  a 
CEP)  . 

4)  Time-late  statistics  (defined  as  the  time  that 
elapses  between  collection  and  processing  of 
reconnaissance  data). 

Values  of  these  performance  parameters  are  computed 
for  each  of  four  data  processing  levels.  The  levels  are: 

1)  Real-time  and  near-real  - time  processing  of 
reconnaissance  data  on  board  the  aircraft. 

2)  Keying  (flagging)  of  data  which  contains  relevant 
target  information  and  transmission  of  this 
keyed  data  to  the  ground  support  station  for 
early  interpretation. 

3)  Keying  of  relevant  data  but  no  transmission  of 
this  data  to  ground  support  station. 

4)  Conventional  processing  (that  is  , no  keying  of 
relevant  data  and  no  transmission  to  ground 
support  station.) 

Thus  the  model  user  can  analyze  increased  system  performance 
as  a function  of  the  type  of  data  processing  used  by  the 
reconnaissance  system. 


m 


mi 


*\ 


1-2 


A 


\ 


% *.  1 


1 i 

■ 


w. 


/v: 


Uv 


•T,"  % V _■"»  ffV  ."V 


1.2  OVERVIEW  OF  IMPROVEMENTS 

Analytics  was  tasked  to  make  improvements  in  the 
AIRS  Model  in  five  major  areas,  as  follows: 


1)  Develop  mathematical  techniques  for  measuring 
the  ability  of  a reconnaissance  system  to 
provide  detectable,  identifiable,  and  localizable 
information  concerning  target  complexes  as 
well  as  individual  targets.  (Before  this  im- 
provement, the. model  considered  only  targets, 
not  target  complexes.) 


2)  Develop  mathematical  techniques  which  can  be 
used  to  measure  the  effectiveness  of  slewing 
optical  sensors  aboard  reconnaissance  aircraft. 

3)  Improve  the  sensor  models  by  inclusion  of 
additional  equipment,  environmental,  and 
operational  parameters. 

4)  Develop  algorithms  that  can  be  used  to  reduce 
the  outputs  of  the  AIRS  performance  simulation 
to  meaningful  measures  of  system  performance. 

5)  Make  such  additions  and  changes  as  are  necessary 
to  provide  additional  AIRS  outputs  which  will 
enable  the  user  to  examine  why  sensors  behave 

as  they  do  on  a micro-level. 
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Task  Cl)  is  the  topic  of  Section  II,  task  (2)  is 
considered  in  Section  III.  Section  IV-  discusses  work  per- 
formed under  task  (4).  Tasks  (3)  and  (5)  are  covered  in 
Sections  V and  VI,  respectively.  Section  VII  provides  the 
AIRS  user  with  instruction  for  the  use  of  the  main  computer 
programs  which  comprise  the  AIRS  Model. 


SECTION  II 


ANALYSIS  OF  TARGET  COMPLEXES 


2 . 1 INTRODUCTION 

The  "Target  Complex"  problem  is  treated  in  re- 
sponse to  a question  highly  relevant  to  air  reconnaissance. 
The  problem  is  to  determine  how  additional  information 
about  specific  targets  can  be  inferred  from  observed  data 
provided  by  the  reconnaissance  system  on  other  targets  in 
the  surrounding  area.  The  approach  pursued  in  this  project 
is  logically  analogous  to  the  thought  processes  of  the  human 
interpreter.  Namely,  if  it  is  concluded  that  each  clement 
of  a group  of  targets  is  functionally  related  as  a member 
of  a target  complex  to  the  other  elements  of  that  group 
(on  the  basis  of  a p r i o r i knowledge  of  the  nature  of  target 
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complexes) , inferences  can  be  made  about  the  identity  of 
more  obscure  targets  which,  by  virtue  of  proximity,  are  also 
members  of  the  complex.  The  following  example  illustrates 
this  more  clearly: 

Suppose  it  is  known  that  a given  type  of  complex  always 
consists  of  targets  arranged  in  the  following  geometric 
pattern  (a  priori  knowledge) : 


Location  1 

Location  2 

Target  Type  1 

Target  Type  1 

Location  3 

Location  4 

Location  5 

Target  Type  3 

Target  Type  2 

Target  Type  3 

Location  6 

Location  7 

Target  Type  1 

Target  Type  1 

Now,  suppose  a photo- interpreter  could  ident i fy  as  well 
as  detect  targets  at  locations  1,2, 3, 5,  and  6,  but 
could  only  detect  (i.e.,  could  not  identify)  targets  at 
locations  4 and  7.  Once  he  realized  he  was  examining 
a known  pattern  (i.e.,  observes  the  complex)  he  can 
make  reasonable  decisions  about  the  types  of  targets 
at  locations  4 and  7 (i.e.,  make  inferences  about  their 
identities) . 


Thus,  the  objectives  of  this  project  arc  necessarily 
two-fold.  The  first  is  to  determine  whether  a target  complex 
is  observed,  and  the  second  is  to  use  this  information  to 
upgrade  the  identif iabi 1 ity  for  each  target -member  of  the 
complex . 
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2.1.1  Definition  of  a Complex 

A target  complex  is  defined  to  consist  of  one  or 
more  types  of  targets;  however,  there  may  be  one  or  more 
targets  within  a type.  Three  kinds  of  complexes  have  been 
studied:  (1)  Anti-aircraft  artillery  (AAA)  sites,  (2) 

Surface-to-Air  Missile  (SAM)  sites,  and  (3)  truck  convoys. 

AAA  and  SAM  sites  are  modeled  as  consisting  of 
up  to  four  target  types  (the  sum  of  the  number  of  targets 
in  the  four  types  is  constrained  to  be  less  than  21).  For 
example,  an  AAA  site  might,  be  modeled  as  consisting  of 
artillery,  radar,  power  vans,  and  computer  vans  --  there 
might  be  4 artillery  emplacements,  2 radars,  2 power  vans, 
and  1 computer  van. 

Convoys  are  modeled  as  consisting  of  only  one 
target  type  --  trucks.  There  may  be  up  to  20  trucks  in  a 
convoy. 

AAA  and  SAM  sites  are  modeled  as  being  "compact" 
complexes.  This  means  that  all  targets  comprising  the  complex 
are  considered  close  enough  in  geographical  proximity  that 
at  any  time  the  probability  that  any  target  is  obscured  by 
terrain  or  clouds  is  the  same  as  the  probability  that  any 
other  target  in  the  complex  is  obscured  by  terrain  or  clouds. 

On  the  other  hand  convoys  are  modeled  as  "dispersed" 
complexes.  This  means  that  the  trucks  which  comprise  the  con- 
voy are  considered  to  be  separated  geographically  to  a degree 
sufficient  to  allow  each  truck  to  be  considered  independent 
of  the  other  trucks  with  respect  to  clouds  and  terrain  blocking. 

2.1.2  Definition  of  Complex  Observability 

A complex  is  considered  observed  when  sufficient 
data  is  developed  by  the  system  on  the  targets  which 
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constitute  the  complex  to  verify  that  a complex  really 
exists.  Since  information  is  supplied  in  a probabilistic 
sense  (i.e.,  detectabilities  and  identif iabilities  of  the 
individual  targets) , a complex  can  be  observed  only  with 
some  probability;  this  probability  is  called  the  observ- 
ability of  the  complex. 


2. 1.2.1  Compact  Complexes.  Surface-to-Air  Missiles  (SAM) 
and  anti-aircraft  artillery  (AAA)  complexes  are  similar  in 
that  they:  (1)  consist  of  several  target  types,  (2)  may  con- 

tain several  targets  of  each  type,  (3)  are  arranged  in 
geographical  areas.  A complex  of  this  category  will  contain 
g (from  one  to  four)  target  types,  and  there  will  be  Nj 


targets  in  type  T^.*  1^  j is  defined  as  the  probability 

(conjoined  with,  not  conditional  on,  detection)  of  identifying 
the  jth  target  of  type  T^,  where  j ranges  from  1 to 
and  h ranges  from  1 to  g. 

A decision  rule,  shown  in  the  following  form,  has 
been  established  by  the  user  for  determining  whether  or  not 
the  complex  is  observed. 


TARGET  TYPE 


Decision 

Function 


1,2  • 


is  constrained  to  be  no  more  than  20  because 
of  practical  limitations. 
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Under  the  decision  rule,  a complex  is  observed  if  any 
decision  function  is  satisfied.  Decision  function  1 is 
satisfied  if  at  least  ^ targets  of  type  and 

Vj  2 targets  of  type  T2  and  . . . and  p targets  of 
ype  T are  identified;  the  other  decisio 


type  T 


on  functions  are 


analogous . 

Let  D^  k be  the  probability  of  identifying  at 
least  k targets  of  type  (k  ranges  from  0 to  N^)  for  a 

given  decision  function.  The  method  for  computing  the 
Dj,  k's  is  given  in  Paragraph  2. 2. 2.1.  The  general  expres- 

_ • r r.  _ _ a • ..  . 


sion 


for  satisfying  decision  function  l is: 


g 

" Dk  h 
h=l  K’n 


(2.1) 


where  k is  The  method  used  for  computing  whether 

the  complex  is  observed,  that  is,  whether  any  decision 
function  is  met,  is  presented  in  Paragraph  2.3. 

2. 1.2. 2 Dispersed  Complexes.  In  a convoy  (the  only  dis- 
persed complex  modeled)  there  is  only  one  target  type, 
namely,  trucks.  Therefore,  the  decision  rule  for  the  complex 
consists  of  a single  decision  function,  and  that  function 
has  only  one  element,  ^ . That  is,  if  at  least  ^ 

trucks  are  identified,  the  complex  is  considered  to  be 


observed.  The  value  of 
parameter . 


is  a user  established 


2 . 2 TREATMENT  OF  COMPACT  COMPLEXES 

The  EXECUTIVE  program  of  AIRS  treats  cadi  target 
independently;  in  particular,  each  target  at  a given  time 
has  a probability  of  being  masked  by  terrain  and  a proba- 
bility of  being  obscured  by  clouds  which  is  independent  of 


any  other  target,  no  matter  how  close  they  may  be  geograph- 
ically. As  already  pointed  out,  compact  complexes  assume 
dependence  of  terrain  and  clouds  for  all  targets  within 
the  complex.  Hence,  a method  was  devised  which  allows 
recomputing  the  identif iabilities  of  each  target  in  the 
complex  assuming  dependence  and  using  as  input  only  the 
independent  identif iabilities  as  computed  by  EXECUTIVE. 
Paragraph  2.2.1  shows  how  this  is  accomplished. 

Once  these  new  identifiabilities  are  computed, 
the  observability  of  the  complex  can  be  computed  by  using 
the  decision  function  for  that  complex.  The  exact  methods 
used  to  compute  the  observability  are  given  in  Paragraph  2.2.2. 

Once  the  observability  of  the  complex  is  knoim, 
the  identifiabilities  of  each  of  the  targets  that  comprise 
the  complex  can  be  upgraded.  This  is  discussed  in  detail  in 
Paragraph  2.2.3. 

2.2.1  Recomputation  of  Target  Identifiabilities 

For  each  target  and  each  processing  level,  EXECUTIVE 
produces  the  probability  of  detection  given  that  there  is  no 
terrain  shadowing,  weather  is  perfect,  and  the  equipment  is 
functioning  properly.  Let  this  probability  be  denoted  by 
BDET^  y where  i is  the  processing  level  and  j is  an 
index  for  all  targets  in  the  complex.  EXECUTIVE  also  com- 
putes a total  probability  of  detection  which  takes  terrain, 
weather,  and  equipment  into  consideration.  This  quantity, 
for  the  jth  targets  in  the  complex,  at  processing  level  i 


f vm 


l^inrvvT.Trvty 


is  a degrading  factor  which  is  the  fraction  of  times  that 
terrain,  weather  and  equipment  parameters  do  not  interfere 
with  detection. 

For  compact  complexes  it  is  assumed  that  there  is 
dependence  among  these  parameters  for  each  target  in  a given 
complex.  Consequently,  an  average  degrading  factor  appli- 
cable over  all  levels  and  targets  is  needed.  This  average, 
denoted  as  DF,  is  found  as  the  average  of  the  individual 


ratios : 


I I 


_ i=l  3=1 
DF  = 


TDET . . 

.».  1 

BDET.  . 


(2.2) 


where  N^,  is  the  total  number  of  targets  in  the  complex. 

The  quantity  DF  is  now  used  to  calculate  the 

modified  total  identif iability  TIDENT.  . for  target  j, 

1 > J 

level  i (i  = 1,  ...»  4 ; j = 1,  ...,  N^)  in  the  following 

manner.  EXECUTIVE  computes  the  conditional  identif iability 
CIDENT.  . , the  probability  of  identification  given 

detection,  not  taking  into  account  terrain,  weather,  and 
equipment  conditions.  Clearly, 

(CIDENTi  . ) (BDET.  ^) 

gives  the  probability  of  identification  conjoined  with 
detection  and  assuming  no  terrain  blocking,  adverse  weather 
conditions  or  equipment  failure.  The  overall  probability  of 
identification,  conjoined  Avith  detection  and  degraded  by  DF 
is  expressed  by: 

TIDENT!  = (CIDENT.  . ) (DF)  (BDET.  .)  (2.3) 

1 > J 1 » 3 1 * 3 


• * ' J 


Zj  ia 


Si^la  rly,  TDIiT.  . is  modified  by  the  following  equation 

.....  1»J  _ . . . . , . 


(the  modified  value  of  TDIiT-  . 

1 9 J 


being  denoted  by  TDET ' - -): 

1 > 3 


TDIiT'.  . = (DF)  (BDET.  .) 

hj  i,3 


(2.4) 


2.2.2  Computation  of  Complex  Observability 

The  first  step  in  computing  the  complex  observ- 
ability is  to  compute,  for  each  target  type,  the  probability 
that  at  least  k targets  of  that  type  arc  identified,  where 
k ranges  from  zero  to  the  number  of  targets  of  that  type. 
This  is  accomplished  in  the  manner  described  in  Paragraph 
2. 2. 2.1.  Paragraph  2. 2. 2. 2 then  gives  the  method  for 
computing  the  probability  that  at  least  one  of  the  decision 
functions  for  the  complex  is  met,  i.e.,  the  observability 
of  the  complex. 


2.22.1  Computation  of  the  Probability  of  Identifying 
nt^gast  k Targets  out  of  a Total  of  for  a Given 

Target  Type.  This  problem  can  be  reworded  as:  Compute 

The  probability  of  realizing  at  least  k out  of  Nj  events 
|iven  that  the  events  do  not  occur  with  the  same  probabilities 

This  cannot  be  computed  in  a straightforward  manner* 
because  the  computations  are  exceedingly  long  and  many 
intermediate  computations  are  required,  thus  the  accuracy 
is  lost  when  they  are  combined. 

In  order  to  accurately  calculate  by  computer  the 


probability  of  realizing  k of 


N 


h 


events,  a short 


Sec  the  discussion  in: 


Feller,  William,  An  Introduction  to  Probability 
Theory  and  Its  Applications  , To  1 ume  T , Tnd~  EcTi  t i on  , 

1 9 b 5 , Wiley  and  Sons,  New  YorF,  New  York,pp.  SS-iUO. 


uhh  wj  u*  VVVV  V '*•  V ■ •-■  •.■  -."  j .■  -.' 


algorithm  was  developed.  This  algorithm  may  be  represented 
by  the  two  general  equations  shown  below.  The  calculations 
arc  performed  iteratively,  where: 

the  iteration  index 

the  number  of  targets  identified 

the  probability  of  occurrence  of 
the  j th  event 


1 

k 

P'V 


a-PCEj ) 


p(Ek,P 


the  probability  that  the  j th  event 
does  not  occur 

the  jth  iteration  of  the  probability 
that  exactly  k events  occur 


N. 


^ = the  maximum  number  of  events  possible. 


p (E0>  j+i}  = P (E0>j)  ^-p  (Ej+l^  (2-5) 

p (Ek>  j+i>  = p (Ej+i)  pCEk.i>j)  + C1-p(Ej+i))p(EkJ) 

(2.6) 

k ranges  from  0 to  j+1 
j ranges  from  0 to  N^-l 


and  P (E^  q)  is  set  to 


1 when  k = 0 
0 when  k > 0 


A formal  proof  of  this  algorithm  is  given  in 
Appendix  A. 

Applying  this  technique  to  each  target  type  in 
the  complex  leads  to  a set  of  probabilities: 

k = 0,1,  ....  N 


k,h  { 


h 


h - 1,2, 


ij 
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where 

Pj,  ^ = the  probability  that  k targets  of 

* type  h are  identified 

= the  number  of  targets  of  type  h 
g = the  number  of  target  types 

Another  useful  statistic  can  now  be  computed  in 
a straightforward  manner;  namely,  the  probability 

that  at  least  k targets  of  type  h are  identified.  It 
can  easily  be  seen  that: 

N. 

“ (k  = 0,1,  N,J 

Dk,h  TT  Pn’h  (h  = 1,2,  ...,  g) 


2. 2. 2. 2 Probability  That  at  Least  One  of  the  General 


Criteria  Sets  is  Satisfied.  Suppose  for  a given  complex, 
the  following  decision  rule  is  given: 


TARGET  TYPE 


Decision 

Function 


2,2  * 


l , 2 


where  ^ is  as  defined  in  Paragraph  2. 1.2.1.  Reitcratins 

decision  function  1 says  that  at  least  V,  . of  target  type 
T ^ and  2 of  target  type  1’2  and  . . . and  of 

target  type  T must  be  detected  for  the  complex  to  be 

o 


2-10 


observed  under  that  decision  function.  The  other  decision 
functions  are  analogous.  The  complex  observability  then 
becomes  the  probability  that  the  complex  is  observed  under 
decision  function  1,  or  decision  function  2,  ...  or 
decision  function  l. 

Suppose  the  probability  of  observing  the  complex 
under  each  rule  is  known,  how  does  one  compute  the  overall 
observability?  The  answer  is  not  clear.  Let  be  the 

probability  that  the  complex  is  observed  under  decision 
function  l.  At  first  glance,  it  would  seem  that  the 
observability,  denoted  by  OB  would  simply  be 

J. 

OB  = 1 - I f(l  - W.) 

i*l  (2.8) 

However,  closer  inspection  reveals  that  this  is  not  the  case, 
the  reason  being  that  the  rules  are  not  independent.  An 
example  might  show  this  dramatically;  consider  a two-target 
type,  two-decision-function  rule: 

TARGET  TYPE 

Decision 

Function  T^  T2 
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the  probability  that  the  least  stringent  decision  function 
is  met. 


Another  approach  would  be  to  let  OB  be  found  as: 
OB  = MAX  [ Wx , W2 , ...  W£  ] 


(2.9) 


Though  this  does  not  have  the  problem  discussed  above,  it 
tends  to  understate  the  observability  by  not  taking  into 
account  the  fact  that  the  complex  could  be  observed  under  a 
non-redundant  decision  function. 


Many  alternatives  were  considered,  but  the  one 
chosen  was  an  iterative  technique  developed  by  Analytics 
which  yields  exact  answers.  It  is  quite  involved  mathema- 
tically, and  if  the  reader  wishes  to  examine  the  mathemat ical 
development  and  proof  of  exactness,  refer  to  Appendix  B. 

It  is  sufficient  here  to  specify  that  the  observability  of 
the  complex  in  level  i will  be  denoted  throughout  the 
remainder  of  this  paragraph  by  OB^. 


2.2.3 


Upgrading  the  Identifiability  of  the  Individual 

Targets 


This  paragraph  illustrates,  once  the  observability 
of  a complex  is  known,  how  the  identif iabilities  of  the 
individual  targets  which  comprise  that  complex  are  upgraded. 
In  order  to  show  these  computations,  several  factors  are 
needed  which  have  already  been  computed,  namely: 


The  total  identifiability,  PI.  .,  of  target 


i,3 


j level  i.  (This  is  just  TIDliN'Ih  ^ as 
computed  in  Paragraph  2.2.1). 


The  total  detectability,  I’D.  .,  of  target  j 

i,3 

level  i.  (This  is  just  TDPT'.  . as  computed 
in  paragraph  2.2.1). 
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The  observability  OB^,  the  complex  at  level  i, 
(Computed  in  Paragraph  2. 2. 2.  2). 


In  addition,  a quantity  called  OB  ^ ^ , must  be  computed.  This 
quantity  is  the  probability  of  observing  the  complex  given 


that  target  j is  not  identified  (i.e.,  given  PI-  . = 0) 

1 > J 

at  level  i.  OB.  . is  found  for  each  target  (j  **  1,  NT) 

at  each  level  (i  = 1,  4)  by  temporarily  letting  PD-  . 

1 y J 

equal  zero  and  recomputing  the  observability  by  the  algo- 
rithm discussed  in  Paragraph  2. 2. 2. 2. 


Suppose  the  complex  has  been  flown  over  (or  near) 

some  number  of  times.  Then  (1  - PI.  .)  represents  the 

1 > J 


percentage  of  times  the  target  was  not  identified  through 
its  own  signature.  This  percentage  should  decrease  by  pre- 
cisely the  number  of  times  the  target  could  be  identified 
through  examination  of  other  targets  in  the  complex  when  the 
target  itself  could  not  be  identified.  In  other  words,  the 
upgraded  identif iability  of  a target  is  the  sum  of  two 
factors:  (1)  the  target's  own  identif iability  without 

regard  to  whether  the  complex  was  observed,  and  (2)  a term 
which  gives  the  identifiability  of  the  target  due  to  the 
complex  observability  when  it  is  not  identif i.able  by  its 
own  signature. 


The  first  factor  above,  for  target  j level  i 


is  PI.  .. 

i.J 


To  develop  the  second  factor,  it  is  assumed  that  a 
target's  identifiability  can  be  upgraded  whenever  the  complex 
has  been  observed  and  the  target  itself  is  detected  with  some 
non-zero  probability. 
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Drawing  these  factors  together  and  expressing  the 

results  mathematically,  PI'.  the  upgraded  identif iability 

1 » J 

of  target  j level  i is: 


PI' . . = PI.  . + (PD.  . - PI.  .)  OB' . . 

i»3  1 i,3  i,3  . 1,3 


(2.10) 


2.3  TREATMENT  OF  DISPERSE  COMPLEXES 

Since  all  targets  within  a disperse  complex  are 
considered  as  independent  with  respect  to  terrain  and 
weather  considerations,  the  identif iabilities  computed  by 
the  EXECUTIVE  program  for  each  truck  in  the  complex  can  be 
used  without  modification  to  ascertain  whether  the  convoy 
(the  only  modeled  disperse  complex)  was  observed.  Paragraph 
2.3.1  discusses  how  the  observability  is  computed. 

Once  the  observability  of  the  complex  is  known, 
the  identifiabili ties  of  the  targets  (i.e.,  trucks)  which 
comprise  the  complex  can  be  upgraded.  The  procedures  for 
accomplishing  this  is  discussed  in  Paragraph  2.3.2. 


2.3.1 


Computation  of  the  Observability  of  the  Complex 


The  EXECUTIVE  program  provides  the  probabilities 

for  (total)  identifying,  by  level,  each  of  the  targets  in 

the  complex.  Let  PI.  . be  the  probability  of  identifying 
. i , J 

the  jth  target  (j  =1,  2,  ...  N^,;  Ny  = the  number  of 

targets)  at  level  i.  Using  the  technique  given  in  Paragraph 

2. 2. 2.1  the  probability  of  detecting  at  least  k targets 

out  of  the  Ny  in  the  complex  is  found  for  level  i;  denote 

this  quantity  by  D^  ^ (k  = 0 , 1 , . . . Ny ; i = 1,  ...,  4). 

Now  the  decision  rule  for  whether  a disperse  complex  is 

observed  consists  of  a single  number,  namely,  the  minimum 

number  of  targets  that  must  be  identified  before  the  complex 

can  he  observed.  Let  this  number  be  k...  . Then,  the 

Min  ’ 
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probability  of  observing  the  complex  at  level  i,  OB.  is 


just  the  probability  that  at  least  k 


identified : 


Min 


targets  are 


OB.  * D. 

1 kMin  * 1 


(i  = 1,  ....  4) 


(2.11) 


2.3.2  Upgrading  the  Identif iabilities  of  the  Targets 
The  upgrading  algorithm  for  disperse  targets 


functions  in  the  same  way  as  that  for  the  compact  complexes 
(see  Paragraph  2.2.3).  Let 


PI-  . = the  total  identif iability  of  target  j 
*•'  level  i * 

PD.  . = the  total  detectability  of  target  j 
level  i * 


OB. 

l 


the  observability  of  the  complex  at 
level  i. 


Using  the  technique  discussed  in  Paragraph  2.2.3, 


OB ' . , the  observability  of  the  complex  given  that  the 


i » j ’ 


identif iability  of  target  j is  zero,  at  level  i,  is  computed. 


Finally  PI'.  . the  upgraded  identifiability 
1 > J 

for  target  j,  level  i,  can  be  written  as: 


PI 


' - - = PI . - + f PD . . - PI . . ) OB' . . 

1 » j 1 » j 1,1  i , r 1,1 


(2.12) 


Since  disperse  complexes  assume  independence,  their 
quantities  arc  not  modified  as  in  the  case  of  compact 
compl exes . 
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2.4  COMPUTER  IMPLEMENTATION 

The  program  which  implements  the  complex  model 
consists  of: 

(1)  A main  program,  COMPLEX,  which  reads  in 
user  parameters  and  data  from  the  AIRS 
EXECUTIVE  routine  and  rearranges  the  data 
for  use  by  the  COMPLEX  subroutines. 

(2)  A subroutine,  SAMAAA,  which  computes  the 
observability  of  compact  complexes. 

(3)  Utility  subroutines  PREPAR  and  OBS  called 
in  by  SAMAAA. 

(4)  A utility  routine,  MULTI,  called  in  by 
PREPAR  and  CONVOY. 

(5)  Utility  function,  PALST,  called  in  by  OBS. 


(6)  A subroutine,  CONVOY,  which  is  used  for  the 
only  disperse  complex,  the  convoy.  CONVOY 
also  utilizes  MULTI. 

Figure  2-1  presents  a general  flow  diagram  of 

COMPLEX. 


2.4.1 


The  Main  Program,  COMPLEX 

The  purpose  of  the  main  program  is  to: 

(1)  Read  the  decision  rules  for  SAM,  AAA, 
and  convoy  complexes. 

(2)  Read  the  target  types  which  comprise  the 
SAM,  AAA,  and  convoy  Complexes. 

(3)  Read,  for  each  complex  to  be  analyzed, 
the  complex  type  number. 

(4)  Read,  for  each  complex,  the  target  number 
of  the  targets  which  comprise  the  complex, 

(5)  For  each  target  in  the  complex,  determine 
by  level: 

(a)  The  basic  detectability 

(b)  The  total  detectability 

(c)  The  conditional  identif iability 

(d)  The  total  identif iability 


The  main  program  develops  these  statistics  before 
(1)  they  are  modified  to  take  into  account  correlation  of 
terrain  and  cloud  cover  for  targets  in  compact  complexes, 
and  (2)  before  they  are  upgraded  as  a function  of  the 
observability  of  the  complex  of  which  the  targets  are 
members . 

Decision  functions  arc  read-in  via  card  input. 

The  following  illustrates  a decision  rule  consisting  of  two 
decision  functions  for  three  target  types. 


DECISION 

FUNCTION 


TARGET  TYPE 


2-17 


A 


' v*  o *.*  ^ • * 


Decision  function  1 requires  that  at  least  0 of  type  23  and 
at  least  3 of  type  4 and  at  least  2 of  type  12  must  be  de- 
tected for  observation  to  occur.  Decision  function  2 is 
similarly  interpreted.  For  overall  complex  observability, 
decision  function  1 or  decision  function  2 must  be  met. 

In  the  case  of  a convoy  the  decision  rule  is 
simple  --a  single  number.  The  probability  that  this  number 
(or  more)  of  trucks  are  identified  is  the  probability  that 
the  convoy  is  observed. 

The  target  types  which  comprise  a given  complex 
arc  found  from  the  SCENARIO  input  deck.  For  example,  in  the 
SCENARIO  input,  target  type  19  might  be  defined  as  power  vans, 
target  type  1 as  SAM  launchers,  etc. 

The  complex  type  and  the  targets  which  make  up  the 
complex  are  also  user  inputs.  Complex  type  1 is  defined  to 
be  a SAM  site;  type  2 an  AAA  site,  and  type  3 a convoy. 

Target  numbers  of  the  targets  which  comprise  a complex  are 
also  read  in  from  cards.  The  way  the  user  determines  the 
target  number  is  to  count  the  position  of  the  target  card 
(for  the  relevant  target)  in  the  SCENARIO  input  deck;  thus, 
the  first  target  card  is  for  target  number  1,  the  second  for 
target  number  2,  and  so  on. 

For  each  target  in  the  complex,  the  main  program 
reads  the  tape  (Output  Tape  #1)  prepared  by  EXECUTIVE.  Since 
the  main  concern  is  in  how  well  the  overall  aircraft  system 
did  for  a given  target,  all  records  on  the  tape  relating 
to  how  a particular  sensor  performed  against  that  target  are 
skipped . 

The  EXECUTIVE  tape  is  organized  on  a pass/target 
basis;  therefore,  the  main  program  must  aggregate  identifi- 
abilitics  and  detectabilities  over  all  passes.  Consider  the 


basic  detectability  at  level  i of  target  j,  pas  k;  call 
this  B.  . , ; further,  let  B'.  . represent  the  probability 
of  (basic)  detecting  target  j over  all  passes,  then: 


...-iV,, 

p=l  1,J’] 


where  there  are  pT  passes.  This  reflects  the  fact  that  for 
a target  to  be  detected,  it  must  be  detected  only  on  1 of  the 
P-p  passes.  The  total  detectability,  conditional  identifi- 
ability,  and  total  identifiability  are  computed  in  an 
analogous  manner. 

Figure  2-2  is  a detailed  flowchart  of  the  function 
of  the  main  program,  COMPLEX. 
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as  they  pertain  to  particular 
sensors) 


Initialize  Level  (13) 
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ditional and  total  identi- 
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Figure  2-2 
(Conti nued) 
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For  each  target,  form: 

(1)  Pr  (conditional  identification) 

(2)  Pr  (total  identification) 

These  are  just  1 - (quantities 

computed  in  box  (16)  for  respective 
target) 


Is  the  complex  type  a compact 
one  (i.e.,  SAM  or  AAA)? 


Divide  total  identif iability  by 
conditional  identif iabil ity  to 
yield  a terrain/weather/equipment 
factor  for  each  target 


Call  Convoy) ( 


(Call  SAMAAA) 
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Figure  2-2 
(Continued) 
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SAMAAA 


The  SAMAAA  subroutine  is  called  in  by  COMPLEX 
whenever  the  complex  being  analyzed  is  either  a SAM  or  AAA 
(Complex  Type  1 or  2).  The  complex  type  is  used  as  an 
argument  of  the  subroutine  while  all  other  inputs  are  through 
common  statements. 


SAMAAA  performs  three  functions: 


CD 


(2) 


An  average  terrain- shadowing , equipment -up , 
and  weather-degradation  factor  for  all 
levels  is  computed  and  used  to  find  the 
total  identif iability  and  the  total  dectect- 
ability  for  each  target  in  the  complex. 

The  probability  of  observing  the  complex 
is  computed  for  each  level  and  the  results 
are  printed. 

The  targets  are  grouped  by  type  and,  for 
each  level,  the  original  probability  of 
identification  is  printed  along  with  the 
augmented  probability  of  identification. 

SAMAAA  calls  in  PREPAR  and  OBS  for  functions  (2) 
and  (3),  respectively,  to  calculate  the  probability  of 
observing  the  complex  at  each  level.  A flowchart  of  the 
SAMAAA  routine  is  shown  in  Figure  2-3. 


(3) 


2.4. 2.1  PREPAR.  This  subroutine  prepares  data  for  the 


OBS,  which  implements  the  decision  function  algorithm 
(Paragraph  2. 2. 2. 2)  to  determine  the  probability  of  ob- 
serving the  complex  at  each  level.  PREPAR  lias  no  arguments; 
all  input  and  output  is  through  common. 
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For  each  level,  PREPAR  uses  the  total  identi- 
fiability  of  each  target  and  computes  the  probability  of 
identifying  at  least  k (k  = 0,  1,  2,  NjJ  targets 

of  each  type. 

PREPAR  calls  MULTI  once  for  each  level  for  each 
target  type.  A flowchart  of  PREPAR  is  given  in  Figure  2-4. 

2. 4. 2. 2 MULTI . This  subroutine  has  as  inputs  N^,  the 

total  number  of  targets  of  a given  type  in  the  Complex,  and 
PD,  an  array  that  comprises  the  probabilities  of  identifi- 
cation of  each  of  the  targets.  The  probability  of 

identifying  at  least  k »(k  = 0,  1,  ...,  N^)  of  the 
targets  is  generated  as  an  array  called  PATLST.  The 
computation  is  based  on  the  iterative  technique  described 
in  Paragraph  2. 2. 2.1.  Figure  2-5  is  a flowchart  of  MULTI. 

2. 4. 2. 3 OBS . Subroutine  OBS  implements  the  decision 
function  algorithm  (Paragraph  2. 2. 2. 2).  OBS  uses  the 
results  of  PREPAR  to  calculate  the  probability  of  observing 
the  complex.  The  type  of  complex,  T^ , is  the  input. 
PROBS(i)  is  an  output  array  of  the  probability  of  com- 
plex observation  at  each  level,  i.  OBS  utilizes  the 
function  PALST  to  compute  PROBS(i).  A flowchart  of  OBS 

is  given  in  Figure  2-6. 
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provided  the  total  identifiability  and  total  detectability 
at  each  processing  level  for  each  truck  in  the  convoy. 

For  each  level,  CONVOY  calls  upon  MULTI  (Paragraph 
2. 4. 2. 2)  to  compute  the  probability  that  at  least  k trucks 
are  identified,  where  k is  the  minimum  number  of  trucks, 
^Min’  that  roust  be  identified  before  the  complex  can  be 
considered  as  observed  (k^.^  is  specified  by  the  decision 
rule  V^  for  convoys).  Thus,  the  observability  of  the 
complex  is  found  for  each  level. 

Next,  again  on  a processing  level  basis,  CONVOY 
considers  each  truck  in  turn.  It  sets,  temporarily,  the 
identifiability  of  the  truck  being  considered  to  zero  and, 
via  MULTI,  computes  the  probability  that  the  complex  is 
observed  given  that  the  identifiability  of  the  truck  is 
zero.  Using  this  probability,  the  upgraded  identifiability 
for  the  truck  is  computed. 

A flowchart  of  the  functions  performed  by  CONVOY 
in  analyzing  convoy  complexes  is  presented  in  Figure  2-7. 

2.5  INPUTS 

The  following  paragraphs  discuss  the  form  of  the 
card  and  tape  input  to  the  program. 

2.5.1  Card  Inputs 

Table  2-1  gives  the  form  of  the  card  input  deck. 
Cards  1 to  4 define  a decision  matrix  for  a SAM  complex 
and  cards  5 to  8 define  an  AAA  complex  decision  matrix.  The 
following  matrix  is  an  example  of  a decision  rule  for  a SAM 
Complex;  AAA  is  analogous. 
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CALL  MULTI 


Start 


Miscellaneous  (1) 
Initialization 


Compute 
Pr (convoy 
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(4) 



LET 

k - 1 

m 


fed 


(•  i = i + 1 


Temporarily  set  the  total 
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Compute  P'obs  the  Pr(obscrve  (7) 
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CALL  MULTI 


Update  the  total  identifi-  (8) 

ability  for  the  kth  truck,  via 

PI  = PI  + (PD  " PI}  P’obs 
Where : 

P = Pr (Detect,  Total,  the  kth 
truck) 
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Set  Pr(total  identif lability  (9) 
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TABLE  2-1 


INPUT  CARD  DUCK 


Definition 


Forma: 


The  number  of  target  types 
in  a SAM  complex  (maximum  of 
four) . 


For  a SAM  complex,  the  target 
type  number  that  appears  is 
the  Ith  column  of  the  SAM 
decision  matrix. 


The  number  of  logically  "OR" 
connected  decision  functions 
for  observing  a SAM  complex, 
i.e.,  the  number  of  rows  in  SAM 
decision  matrix  (maximum  of  four) . 


For  a SAM  site,  decision  function 
1,  the  min.  number  which  must  be 
identified  of  target  type  in  col. 
I of  SAM  decision  matrix. 


For  a SAM  site,  decision  function 
2,  the  min.  number  which  must  be 
identified  of  target  type  in  col. 
I of  SAM  decision  matrix. 


A T l 


The  number  of  target  types  in  an 
AAA  complex  (maximum  of  four) . 


For  an  AAA  complex,  the  target 


i 414 


type  number  that  appears  in  the 
I tn  column  of  the  AAA  decision 


matrix . 


NRULES  now  becomes  the  number  of 
logically  OR-connected  decision 
rules.  For  observing  an  AAA 
complex,  i.e.,  the  number  of  rows 
in  the  AAA  decision  matrix  (max. 
of  four) . 


R NF 


Table  2-1  - continued 


8 a | DECMAT  (2 , I,  1) 


8b  DECMAT  (2,1,1) 


9 IDECON 


10  NOCPLX 


11  1 ICPXTP 


NOTARG 


12  ' ITAR(I) 

1=1,  NOTARG 


For  an  AAA  site,  decision  function 
1,  the  min.  number  which  must  be 
identified  of  the  target  type  in 
column  I of  AAA  decision  matrix. 


14 


For  an  AAA  site,  decision  function  414 
2,  min.  number  which  must  be  iden- 
tified of  the  target  type  in  columr 
1 of  AAA  decision  matrix. 


The  min.  number  of  trucks  that 
must  be  identified  for  the  convoy 
to  be  observed. 


The  number  of  complexes  which 
are  to  be  analyzed  on  this  run. 
(No  limit  on  this  number.) 


The  complex  type 
1 = SAM,  2 = AAA,  3 = convoy 


The  number  of  targets  (not  types 
in  complex  (maximum  of  ?0  targets 
spread  over  up  to  4 types)* 


The  target  numbers  of  those 
targets  in  the  complex. 


2014 
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Targets  should  be  drawn  only  from  the  type  specified 
in  the  element  decision  matrix.  Targets  of  inappropriate 
type  will  not  affect  observation  of  the  complex,  but 


warranted . 


Complex . 
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each 
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V v V , 


Decision 

Function 


23 


TARGET  TYPE 


9 


4 12 

1 0 3 2 1 

2 5 10  0 

3 2 2 2 2 

MAXTAR(l)  = 4,  i.e.,  there  are  four  target  types.  TARTYP(1,1) 
= 23,  i.e.,  target  type  23  appears  in  the  first  column  of  th 
matrix;  likewise  TARTYP(1,2)  = 4;  TARTYP(1,3)  = 12; 

TARTYP(1,4)  = 9.  Decis ion' function  1 says  that  in  order  for 
the  complex  to  be  observed,  at  least  0 of  type  23  and  3 
of  type  4 and  2 of  type  12  and  1 of  type  9 must  be  identified; 
decision  functions  2 and  3 are  analogous.  In  this  example 
for  the  DECMAT  array,  the  first  subscript  is  1 if  considering 
a SAM  complex  and  2 if  an  AAA  complex;  the  second  subscript 
is  the  column  number,  while  the  third  index  is  the  row  (i.e., 
decision  rule  number)  number.  Note  that  for  a complex  to  be 
observed,  Rule  1,  Rule  2,  or  Rule  3 must  hold,  not  some 
combination  of  them. 


2.5.2  Tape  Inputs 


The  tape  input  is  on  Output  Tape  #1  of  the  EXECUTIVE 
routine.  This  tape  is  discussed  in  detail  in  Paragraph  3.6.1 
of  the  cited  AIRS  report. 

2.6  OUTPUTS 

The  first  page  of  output  presents  the  decision 
matrix  (i.e.,  the  decision  rule)  for  SAM  complexes;  page  2 
is  the  AAA  decision  matrix.  Similarly,  page  3 presents 
the  decision  rule  for  observing  convoys.  Sample  outputs 
are  given  in  Figure  2-S. 
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ON  FIRST  PAGE: 
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DECISION  MATRIX  FOR  SAM  COMPLEX 

TARGET  TYPE 

DECISION  RULE  * 22  21  20 

******************************************** 

1 * 4 1 1 

2 * 2 2 2 

3 * 5 0 0 


ON  SECOND  PAGE: 


DECISION  MATRIX  FOR  AAA  COMPLEX 


TARGET  TYPE 

DECISION  RULE  * 3 6 7 


* * * * * * * * * * * * * * * * * * * * * 


* * * * * * * * * * * * * * ******** 


1 *201 

2 *210 


ON  THIRD  PAGE: 
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FOR  CONVOYS,  AT  LEAST  3 TRUCKS  MUST  BE 
IDENTIFIED  FOR  THE  COMPLEX  TO  BE  OBSERVED 
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Figure  2-8 

Sample  Outputs  --  Decision  Rules 
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For  each  complex  analyzed,  header  information  is 


generated  consisting  of  the  complex  number  (the  first  complex 
analyzed  is  complex  «fl,  the  second  is  complex  #2,  etc.)  and 
the  type  of  the  complex.  Then , for  each  of  the  four  pro- 
cessing levels  the  probability  of  observing  the  complex  is 

* given,  appropriately  labeled.  Within  each  level  individual 
targets  are  grouped  by  type  (in  the  case  of  SAM  and  AAA 

} complexes  --  this  is  not  done  for  convoys,  which  consists 

* 

| of  only  one  target  type).  The  identifiability  of  each 

• target  before  consideration  of  whether  a complex  is  observed 

is  given,  as  in  the  new  target  identifiability  computed 
;•  once  the  probability  of  observing  the  complex  of  which  the 

target  is  a member  is  known.  Sample  outputs  for  SAM,  AAA, 
and  Convoy  Complexes  are  given  in  Figure  2-9. 
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SECTION  III.  PHOTO- SLEWING 


3.1  INTRODUCTION 

The  term  photo-slewing  indicates  the  process 
of  directing  the  orientation  of  reconnaissance  photo- 
graphic devices  with  respect  to  specific  geographical 
points  of  interest.  The  original  AIRS*  model  conceptual- 
izes fixed  photo-sensors,  i.e.,  the  orientation  of  the 
projected  f ield-of -view  with  respect  to  the  aircraft  is 
fixed,  for  a given  mission,  according  to  the  designation 
of  the  user. 


Reported  in  Volumes  I,  IT,  and  III  of  "SIMULATION  OF 
ADVANCED  INTEGRATED  RECONNAISSANCE  SYSTEMS,"  submitted 
to  Systems  Analysis  and  Engineering  Department,  U.S. 
Naval  Air  Development  Center,  Johnsville,  Warminster, 
Pennsylvania  by  Analytics,  Inc.,  under  Contract  No. 
N62269-6S-C-0441 ; April  4,  1969. 
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In  future  reconnaissance  systems  slewing 
capabilities  will  become  increasingly  important,  as  more 
targets  can  be  detected  by  a slewable  camera  than  by  an 
otherwise  identical  fixed-position  camera.  Analytics  was 
tasked  to  improve  the  AIRS  models  to  allow  slewable  cameras 
to  be  modeled. 

The  modeling  technique  chosen  to  implement  this 
slewing  capability  was  to  retain  the  concept  of  a camera 
fixed  in  orientation.  By  modifying  only  relevant  film 
dimensions  and,  in  the  case  of  the  side-looking  frame 
camera,  the  depression  angle  of  the  center  beam,  the  modi- 
fied camera  can  be  made  to  "see"  the  entire  area  over  which 
its  real-life  counterpart  could  be  pointed.  The  field  of 
view  of  this  modified  camera  is  cal  1 ed  the  "expanded  field 
of  view"; the  field  of  view  of  the  original,  non-modificd 
camera  is  referred  to  as  the  "rest  field  of  view." 

Detectability  and  idcnt.i  f iability  are  then 
computed  using  the  expanded  field  of  view,  without  regard 
to  whether  or  not  the  camera  is  actually  slewed.  These 
detectabilities  and  identif iabilities  arc  then  modified 
according  to  the  following  algorithms: 

(1)  If  the  target  would  have  passed  through 
the  field  of  view  of  the  camera  in  its 
unslewed  mode  (a  deterministic  computation), 
slewing  gains  nothing,  and  the  detectabil- 
ities and  identif iabil itics  remain  as 
computed . 

(2)  If  the  target  would  not  have  passed  through 
the  field  of  view  of  the  camera  in  its  un- 
slewed mode,  the  detectability  (identi- 
fiability)  must  be  modified  by  multiplying 
the  computed  detectability  ( ident if iabi 1 i ty) 
by  the  probability  that  slewing  occurred 
for  the  target. 

The  forward-looking  and  side-looking  frame  cameras 
are  modeled  as  being  capable  of  direction  by  preplanning  or 


by  real-time  signals  from  the  ECM,  MTIFLR,  or  FLIR  sensors. 
The  choice  of  these  three  sensors  as  sources  of  slewing  is 
discussed  in  Paragraph  3.4.3.  The  computer  implementation 
has  been  effected  by  additions  to  the  existing  AIRS 
SCENARIO  and  EXECUTIVE  programs. 

Paragraph  3.2  discusses  those  parts  of  the 
slewing  model  which  are  implemented  in  SCENARIO. 

Paragraph  3.3  discusses  the  changes  effected 
in  the  EXECUTIVE  program. 

3.2  SCENARIO 

Paragraph  5.2.1  presents  a discussion  of  the 

geometry  of  the  field-of -view,  in  slewed  versus  unslewed 
mode,  for  the  forward-looking  frame  camera  and  an  expla- 
nation of  the  means  employed  to  achieve  the  larger  field 
of  view.  Paragraph  3.2.2  provides  a similar  treatment  of 
the  side-looking  frame  camera. 

3.2.1  Forward-Looking  Frame  Camera 

In  the  original  AIRS  model,  forward-looking 

frame  cameras  were  not  modelled  as  being  slewable  by  other 
sensors.  The  model  lias  now  been  modified  to  allow  slewable 
forward-looking  frame  cameras  to  be  considered  within  the 
model.  The  modeling  technique  chosen  was  to  continue  repre- 
senting the  camera  as  fixed  in  orientation,  but  to  enlarge 
the  field  of  view  of  the  camera  to  correspond  to  the  entire 
area  over  which  the  slewed  camera  could  be  directed.  However 
resolution  and  other  performance  characteristics  of  the 
slewed  camera  are  the  same  as  those  of  its  unslewed  counter- 
part . 

The  problem  as  appraised  is  to  find  the  ratio  of 
the  area  of  coverage  in  the  slewable  mode  to  the  non- slewable 
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mode.  Since  the  forward  camera  sweeps  a strip  under  the 
aircraft,  this  is  equivalent  to  finding  the  ratio  of  the 
length*  of  the  swaths  swept  out.  This  ratio  can  then  be 
used  to  find  an  effective  film  length,  i.e.,  the  film  length 
which  would  yield  the  path  the  slewable  camera  could  cover. 

Let  f be  the  focal  length  and  X the  film 
length.  Then  a , the  half  angle  of  the  beam  spread 
parallel  to  the  wing  is  given  by: 

a = arctan  (A/2f)  (3.1) 

Let  d be  half  the  swath  length  of  the  field  swept  out: 

d = X tan  a,  (3.2) 

where  X is  the  slant  range  to  the  ground.  X is  given  by: 

X = H/cos  0.  (3.3) 

Where  H is  the  aircraft  height  above  mean  ground  level  and 
0 is  the  angle  between  the  vertical  and  the  uppermost  ray 
of  the  beam.  Equations  (3.1),  (3.2),  and  (3.3)  arc  taken 
from  Figures  3-la,  3-lb,  and  3-2. 

Let  tp  be  the  half  angle  through  which  the 
camera  can  be  slewed,  as  shown  in  Figure  3- lb.  Referring  to 
this  figure,  the  effective  film  length,  X'  , can  be  found 
in  the  following  manner. 

From  equations  (5.1)  and  (3.2): 

d - XX/2 f (3.4) 

From  equation  (3.4)  and  Fibure  3 -lb: 
d’  = (X  + e)  tan  (a  + \p) 


* Length  is  measured  as  parallel  to  the  a/c  wing  axis. 


(3.5) 


It  can  be  assumed  that  * e <<  X. 


Therefore : 

d/d'  = X/f  tan  (a  + ip) 

By  similar  triangles: 

d'/d  = tan  (tt/2  - 0)/tan  a 
= 2f  tan  (a  + ^)/X 
= cot  3/tan  a 


(3.6 


(3.7 


cot  0 = 2f  tan  (a  + ip)  (tan  a)/X 


(3.8 


Then : 

X'  = 2f/tan  0 
= 2f  cot  3 

= ((2f)2/x)  tan  (a  + tan  a (3.9 

which,  by  inspection,  has  the  correct  unit  of  distance. 
Equation  (3.9)  can  be  simplified  to  yield 

X'  = (2f)2/X  tan  (arctan  ( X / 2 £ ) + ip)  • tan  (arctan  X/2f) 

= 2f  tan  (a  + ip)  (3.10 

Empirical  data  indicates  that  within  a certain 
number  of  degrees,  y,  from  the  horizon  (usually  y is 
about  6°),  optical  sensors  are  not  productive.  Therefore, 
we  shall  constrain  6 to  be  not  greater  than  the  slew  angle 
which  would  allow  the  camera  to  sec  within  y degrees  of 
the  horizon.  This  angle  can  be  found  by  referring  to  Figure 
3-3. 


Exactly,  e-  f - e';  e'  = hX  sin  (j  - a.)/  sin  ( ^ - a) 
For  a typical  case  a = 4^,  i p ~ 60°7  X = 2.5  inches, 
f = 16  inches,  e'  = 1.5  inches,  and  e = 14.5  inches 
compared  to  x which  is  several  miles. 
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As  before,  0 is  the  angle  from  the  vertical  to  the  top  edge 
of  the  beam.  6',  as  defined  in  Figure  3-3  is: 

$'  = tt/2  - y (3.11) 

Note  that : 

tan  5 = D/X  = (H/X) (D/H)  = (cos  6)  D/H  (3.12) 

where 


D/H  = (Z2  - Y2>Vh  - (Z2/H2  + Y2/  \{2)h 

= [tan2  (tt/2  - y)  - tan2  e]  **  (3.13) 


Therefore : * 


6 = arctan  ^ cos  0 (tan2  (tt/2  - Y)  - tan20}?5 

\3.14) 

And  finally: 


ii>'  = (-a)  + min  [a  + 4/,  <$]  , (3.15) 


For  a limit  of  6°  from  the  horizon,  6 becomes 
arctan  [(cos  2 0/ . 0109)  - l]?i 


__a  r ~ •f  “ 


where  is  the  user  input  value  of  the  half  angle  of  the 

slew.  Equation  (5.10)  is  then  rewritten  simply  as: 

X'  = 2f  tan  (a  + ty')  (3.16) 


3. 2. 1.1  Computation  of  the  Target  Offset  Distance.  One 
of  the  statistics  that  must  be  provided  to  the  EXECUTIVE 
program  for  each  look  at  a given  target  is  the  offset  distance 
from  the  target  to  the  heading  axis  at  the  time  of  that  look. 

Let : 


F = offset  distance 

f = heading  angle  measured  clockwise  from 
north  (i  . e . , x axis) 

(Yt,Xt)  = the  coordinates  of  the  target 


(Ya,Xa)  = the  coordinates  of  the  aircraft 


Then,  with  the  aid  of  Figure  3-4: 

Ground  distance  to  target,  GR,  is: 

rR  ‘ UXT  ' XA)2  + (YT  ' YA)2f‘ 

Angle  from  north  to  target,  p,  is: 

p = [it/ 2 - arctan  ( (XT  - XA)/(YT  - YA))] 
Angle  from  heading  to  target  = p - f 


Then : 


Gr  sin  (p  - 6) 


(3.17) 
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X axis 


Y axis 

NOTE:  Dashed  areas  show  the 

large  field  or  view , 
i.e.  the  slewed  field 


Figure  3-4 

Target  Offset  Distance 
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3 . 2 . 1 . 2 Computation  of  the  Average  Length  of  Rest  Field 
of  View.  The  rest  field  of  view  is  a trapezoid;  however, 
for  reasons  which  will  be  discussed  in  Paragraph  3.3.1,  a 
necessary  parameter  is  the  average  width  of  the  rest  field 
of  view,  denoted  by  D. 

The  length  of  the  longer  side  of  the  trapezoid  is: 
D,  = 2X  tan  a (3.18) 

Where  X is  computed  by  equation  (3.3)  as: 

X = H/  cos  6 

Let  £ be  the  angle  between  the  aircraft  vertical  and  the 
center  of  the  focal  beam.  Since  0 equals  £ + u,  from 
Figure  3-lc,  D becomes: 

El  = 2H  tan  (a)  sec  (£)  (3.19) 


I 


3 . 2 . 1 . 3 Input  for  Forward-Looking  Frame  Camera  Slewing. 

One  new  parameter  must  be  provided  the  SCENARIO  program  to 
allow  for  slewing  of  forward-looking  frame  cameras,  namely, 
the  half  angle  of  slewing  capability  (in  degrees).  Paragraph 
7.2.1  gives  the  new  form  of  the  SCENARIO  input  deck,  in- 
cluding this  parameter. 


3 . 2 . .1 . 4 Output  Related  to  Forward-Looking  Frame  Camera 
Slewing.  The  modeling  of  the  slewing  capability  of  the 
forward-looking  frame  camera  produces  no  changes  in  the  form 
of  the  SCENARIO  output,  either  hard  copy  or  tape. 
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3.2.2 


Side-Looking  Frame  Camera 


Side-looking  frame  cameras  were  not  modeled  as 
being  slewable  in  the  original  AIRS  model;  this  capability 
has  now  been  added.  As  in  the  case  of  the  forward-looking 
frame  camera  (Paragraph  3.2.1),  the  modeling  technique  chosen 
is  to  consider  the  slewable  camera  as  a fixed-orientation, 
non-directable  camera  having  a field  of  view  which  encompasses 
the  entire  area  over  which  the  slewable  camera  could  actually 
be  directed.  This  is  accomplished  by  modifying  only  physical 
parameters  (e.g.,  film  lengths)  of  the  camera;  performance 
characteristics  (e.g.,  resolution)  remain  unchanged. 

As  in  the  case  of  the  forward-looking  camera,  it 
is  also  assumed  here  that  when  all  other  parameters  are  held 
constant  the  length  of  the  field  of  view  is  proportional  to 
film  length.  The  augmented  film  length,  denoted  by  A'  (see 
Figure  3-5)  is  computed  by  equations  (3.22)  and  (3.23)  where; 


X 

f 

T 


MIN 


MAX 


the  original  film  length 
the  focal  length 

the  angle  between  the  focal  axis  and  the 
aircraft  vertical.  When  the  user  wishes  the 
camera  to  point  to  the  right,  relative  to  the 
direction  in  which  the  aircraft  is  pointed, 

T is  positive.  When  the  camera  should  be 
pointed  to  the  left  (when  looking  along  the 
flight  path)  T is  input  as  a negative  value 

the  minimum  depression  angle  of  slewing 

the  maximum  depression  angle  of  slewing. 

0MIN  and  8^^  are  always  positive.  The 

angles  arc  measured  relative  to  the  horizon 
of  interest  depending  on  whether  the  camera 
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is  right-mounted  or  left-mounted  on  the 
aircraft. . 

Left  Mounted 
(positive) 

| " " a/c  heading 


(positive^  1 Nadir 
_Right  Mounted 

El  = the  signed  magnitude  of  the  quantity 
11/2  ‘ °MIN 

E2  = the  signed  magnitude  of  the  quantity 
7r/2  ' 6 MAX 

a = beam  spread  = arctan  (,\/2f) 

In  order  for  the  upper  edge  of  the  expanded  focal 
beam  to  remain  6°  below  the  horizon,  El  must  not  exceed 
(84°  - a) . For  the  reverse  edge  of  the  expanded  focal  beam 
to  remain  6°  below  the  horizon,  E2  must  not  be  more 
negative  than  (-84°  + a).  V.'hen  El  and  E2  are  set  equal 
to  either  their  original  values  or,  where  they  exceed  the 
computed  limits,  the  focal  axis  is  repositioned  to  bisect 
El  + E2  by  equation  (3.20):  this  repositioned  focal  axis 

is  at  an  angle  denoted  by  T'  from  the  aircraft  vertical. 
Hence , 


T'  = 0 . 5 (El  + E2 ) 


(3.20) 


From  Figure  3-5  it  can  be  seen  that  the  portion 
of  augmented  beam  angle,  a',  above  the  aircraft  vertical  will 
be : 

o'  = T'  + a + 0.5 (El  - E2) 
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a"  = -T'  + a + 0 . 5 (El  - E2) 

and  the  total  beam  angle  is  a'",  given  by: 

a'"  = a'  + a"  = 2a  + El  - E2  i 

and,  since  X'  is  perpendicular  to  the  bisector  of  the 
total  beam  angle  (see  Figure  3-5): 


X'  = 2f  (tan  0.5a'")  (3. 

= 2 f (tan  (a  + 0.5(E1  - E2))j  (3. 

Computation  of  the  Tareet  Offset  Distance.  The 


3. 2. 2.1 


target  offset  distance  is  computed  using  equation  (3.17) 
in  Paragraph  3. 2. 1.1. 

3 . 2 . 2 . 2 Computations  of  the  Right  and  Left  Edges  of  the 
Rest  Field  of  View.  As  will  be  discussed  in  Paragraph 
3.3.1,  the  EXECUTIVE  program  will  need  to  know  the  distance 
from  the  heading  axis  to  the  right-most  edge  of  the  rest 
field  of  view  and  the  similar  distance  to  the  left-most 
edge  of  the  rest  field  of  view.  This  paragraph  will  show 
how  these  distances  arc  computed. 

Whether  or  not  the  camera  is  slewed,  the  actual 
rest  field  of  view  is  delimited  by  T (the  original  angle 
between  the  aircraft  vertical  and  the  focal  axis)  and  a 
(the  half-angle  of  beam  width  corresponding  to  the  original, 
unmodified  film  length).  For  either  right-  or  left-mounted 
cameras,  the  distance  fron  the  heading  axis  to  the  right- 
most edge  of  the  rest  field  of  view  is: 


S , = 11  tan  (T  + a) 


(3.25) 


"▼]  - 


The  corresponding  distance  to  the  left-most  edge  is: 

S,  = II  tan(T  - a) 


3. 2. 2 . 3 


Input  for  Side-Looking  Frame  Camera  Slewing  Model 


Two  additional  parameters  must  be  input  into  the  SCENARIO 
program;  they  are: 

(!)  eMIN  (in  d egrees)  ; defined  in  Paragraph  3.2.2 

(2)  degrees);  defined  in  Paragraph  3.2.2 

New  input  to  SCENARIO  is  described  in  detail  in  Paragraph 
7.2.1  of  this  report. 

3 . 2 . 2 . 4 Output  Related  to  Side-Looking  Frame  Camera 
Slewing . The  modeling  of  the  slewing  capability  of  the 
side-looking  frame  cameras  produces  no  changes  in  the  form 
of  the  SCENARIO  output,  either  hard  copy  or  tape. 


3.3.1 


EXECUTIVE 
Field  of  View 


The  preceding  paragraph  discussed  how  the 
slewing  problem  for  forward  and  side  oblique  frame  cameras 
was  implemented  in  the  SCENARIO  program.  This  paragraph 
discusses  how  the  outputs  provided  by  the  SCENARIO  program 
are  used  in  the  EXECUTIVE  program  to  compute  target  de- 
tectabilities, identifiabiliti cs , and  locali zabilities 
when  slewing  capabilities  are  present. 

Below  are  definitions  of  terms  that  will  be  used: 
A.  = the  distance  of  the  target  from  the  head- 

1 til 

ing  axis  at  the  time  of  the  l look.  It 
is  taken  as  negative  if  the  target  is  to 
the  left  of  the  axis  and  positive  if  to 
the  right . 


A 
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B = the  average  distance  from  the  heading  axis 
to  the  right-most  edge  of  the  rest  field 
of  view.  B is  positive  if  the  right  edge 
of  the  field  of  view  is  to  the  right  of 
the  heading  axis. 

C = the  average  distance  from  the  heading 
axis  to  the  left-most  edge  of  the  rest 
field  of  view.  C is  negative  if  the 
left  edge  is  to  the  left  of  the  heading 
axis . 

Note  that  for  a forward-looking  frame  camera, 

B = D/2  = -C  (Paragraph  3. 2. 1.2); 


for  a side-looking  frame  camera, 

B = and  C = (Paragraph  3. 2. 2.  2). 

The  terms  are  pictorially  illustrated  in  Figure  3-oa  and  3-6b. 
For  each  look: 

if  C < A < B 
otherwise 


Lk  = S 
V 


Aggregating  over  x 


looks  at  a given  target: 


E 


Lk 


h 


If  E _>  0.5,  the  majority  of  looks  at  the  target  fell  with- 
in the  rest  field  of  view.  If  E < 0.5,  the  majority  of 
looks  fell  outside  the  rest  field  of  view,  but  within  the 
expanded  field  of  view.  The  look  would  not  be  recorded  if 
the  target  were  outside  the  expanded  field  of  view. 

Figure  3-7  will  illustrate  different  values  of  E.  It  is 
assumed  that  if  E < 0.5  the  target  cannot  be  detected 


on  any  look  unless  sleAving  occurs;  likewise,  if  E ^ 0 . 5 
it  is  assumed  that  the  target  can  be  detected  on  any  or 
all  looks  even  if  slewing  does  not  occur. 

The  assumption  can  be  justified  by  noting  that 
E is  usually  1 or  0 and  that  only  in  such  very  special 
cases  as  illustrated  in  Figure  3-7  will  E be  between 
these  values. 

3.3.2  Computation  of  Detectability  and  Ident if iabi 1 i ty 

Detectability,  identifiability  and  localizability 
are  initially  computed  assuming  all  looks  are  valid,  i.e., 
using  the  expanded  field  of  view  but  without  regard  to  whether 
slewing  actually  occurred.  These  detectabilities  and 
identifiabilities  may  then  be  modified  to  take  into  account 
the  probabilities  that  slewing  did  or  did  not  occur. 

If  the  camera  would  have  seen  the  target  whether 
or  not  slewing  occurred  (i.e.,  E ^ 0.5),  then  slewing 
provides  no  additional  improvement.  Therefore,  the  detect- 
ability and  identifiability  for  the  target  remains  unchanged. 

If  hoAvever,  the  camera  cannot  see  the  target 
unless  slewing  occurs,  we  have: 

PD ' . = PDt  • F 

where  PD'^  is  the  modified  detectability  for  level  i; 

PD^  the  initially  computed  detectability  for  level  j , and 
Ps  is  the  probability  that  the  camera  was  slewed. 

Likewise 

PI ' . = PI . • P 
1 is 

where  PI ' ^ is  the  modified  identifiability  at  level  i 
and  PI  the  initially  computed  identifiability  at  level  i. 
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Case  2 
E < .5 


Target  Path 


Figure  3-7 

Illustration  of 
Values  of  E Factor 
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I ho  drawings  arc  shown  for  a forward  oblique  frame  corpora; 
the  side  oblique  frame  camera  in  completely  analogous. 
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3.3.3 


Slewing  can  occur  in  one  of  two  ways:  the  target 

can  be  preplanned,  or  the  target  can  be  sensed  by  another 
sensor  which  can  direct  the  slewable  cameras  toward  the 
target.  In  the  preplanned  mode,  the  system  receives  co- 
ordinates of  known  or  suspected  targets  at  the  time  of 
mission  initiation.  For  the  second  mode  the  sensors  must 
be  considered  which  may  reasonably  be  expected  to  pick  up 
a target  in  time  to  direct  the  camera  toward  the  target. 

The  target  must  be  sensed  before  it  comes  abreast  of  the 
aircraft,  for  the  side  oblique  camera  or  while  it  is  still 
well  in  front  of  the  aircraft,  for  the  forward  oblique  camera. 
Further,  the  target  must  be  sensed  in  real  time,  i.e., 

Level  1 of  processing.  The  ECM,  MTIFLR,  and  FLIR  sensors 
meet  these  criteria.  Each  can  detect  targets  in  front  of  the 
aircraft  and  each  operates  in  Level  1.  Because  (see  Para- 
graph 3.2. 8.1  of  AIRS  Volume  I)  FLIR  itself  is  slewable  by 
the  ECM  and  MTIFLR  sensors,  if  FLIR  has  been  slewed,  the 
cameras  also  have  been  slewed,  and  FLIR  can  add  nothing. 

The  slewing  factor  Ps  is  then: 


Where 


Ps  = 1 - R (l-PDx)  (1-PD  ) (l-PDz(l-P's)) 


R 


PD. 


PD. 


y 


PD 


Probability  that  the 
planned . 

Probability  that  the 
target  in  Level  1. 

Probability  that  the 
target  in  Level  1. 

Probability  that  the 
target  in  Level  1. 

Probability  that  the 


target  was  not  pre 
ECM  detects  the 
MTIFLR  detects  the 
FLIR  detects  the 
FLIR  is  slewed. 


and,  from  equation  (3.150)  of  AIRS,  Volume  I. 

P's  = 1 - R (l'PDx)  (l-PDy) 

No  change  need  be  made  in  the  CEP,  because  it  is  not  con- 
ditioned on  detection  or  slewing. 

3-3.4  EXECUTIVE  Program  Input  and  Output 

There  is  no  change  in  the  form  of  EXECUTIVE 
input*  or  output. 


Certain  parameters  are  redundantly  input  to  both  SCENARIO 
and  EXECUTIVE.  The  user  should  note  the  duplication  of 
focal  length  on  input  card  16b  (p . 2-34)  and  input  card 
20  (p.  3-90,  Vol.  I.  Simula  t i on  of  Advanced  Integrated 
Reconnaissance  Systems  (AIRS).) 


SECTION  IV 
EV ALU AT I ON  ROUT I N H 

4.1  INTRODUCTION 

The  Evaluation  (EVAL)  routine  is  a CDC  32/3300 
FORTRAN  program  that  prints  usable  measures  of  AIRS  effect- 
iveness in  a highly  readable  format  by  transforming  data 
contained  on  the  Target/Sensor  Output  Tape  (Executive  Output 
Tape  If  1;  see  Paragraph  3.6.1  of  AIRS,  Vo]  . I and  naragranh 

7. 3. 2.1  of  this  report).  Through  the  use  of  ontions  the 
user  of  EVAL  can  to  some  degree  control  the  type  of  measures 
computed  as  well  as  the  form  of  the  output  printed;  it  is 
through  such  control  that  EVAL  gains  its  flexibility.  Further 
the  user  can  process  the  same  data  tape  repeatedly  using 

different  options  each  time. 
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The  basic  purpose  of  EVAL  is  to  integrate  the 
effectiveness  measures  of  detectability,  identifiability , 
and  localizabili ty  computed  by  EXECUTIVE  over  accumulations 
of  targets  known  as  types , and  further  to  integrate  over 
collections  of  types  known  as  groups . These  measures  are 
integrated  by  sighting  unit , i.e.,  with  each  sighting  of  a 
target  considered  a separate  entity,  and  also  by  target  unit, 
i.e.,  with  the  measures  for  all  sightings  of  a particular 
target  conjoined  into  single  target  measures  before  any 
integration  takes  place.*  Both  target  unit  and  sighting  unit 
data  are  useful  in  evaluating  the  effectiveness  of  AIRS 
systems,  although  sighting  unit  measures  are  especially 
meaningful  when  considering  individual  sensor  performance, 
as  are  target  unit  measures  when  examining  the  system  as  a 
whole . ** 

Considering  each  target  type  in  turn,  the  program 
outputs  effectiveness  measures  accumulated  by  type  (type 
measures)  in  both  target  and  sighting  units.  These  measures 
are  given  for  each  sensor  and  also  combined  over  sensors . In 


* In  a multiple  pass  mission,  the  reconnaissance  system  is 
exposed  to  each  physical  target  more  than  once.  Each  such 
exposure  we  call  a target  sighting  or,  simply  a sighting. 

**  This  is  actually  a deep  Point.  When  target  unit  measures 
arc  considered,  the  system  performance  is  not  adversely 
affected  if  a target  is  poorly  detected  (or  identified,  etc.) 
on  one  leg  of  the  mission  as  long  as  it  has  been  well 
detected  on  another.  This  is  as  it  should  be,  after  all,  the 
purpose  of  a reconnaissance  flight  is  to  detect  targets , and 
less  concern  should  be  given  to  performance  against  individual 
sightings.  On  the  other  hand,  when  the  performance  of  separate 
sensors  is  being  examined,  sensor's  effectiveness  measure 
is  based  on  how  well  it  does  at  every  opportunity;  this  i_s 
affected  if  a sensor  docs  poorly  on  one  sighting  of  a target 
even  if  it  does  well  on  another. 
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addition,  some  break  out  by  processing  level  is  given.  A 
complete  integration  over  all  targets  (both  by  sighting  and 
target  units)  is  then  presented.  Next,  data  for  targets  of 
the  appropriate  types  are  accumulated  into  target  group  data; 
however,  far  fewer  group  measures  are  given  - only  summary 
data  are  tabulated  - and  there  is  no  break  out  by  sensors. 
Lastly,  EVAL  prints  the  beta  measure  of  sensor  redundancy 
both  by  target  and  sighting  unit  for  all  sensores,  integrated 
over  all  targets.* 

The  EVAL  routine  as  discussed  above  has  been  im- 
plemented and  thereby  replaces  and  obsoletes  an  earlier,  more 
limited  version  of  EVAL  (see  Section  IV,  AIRS,  Vol  I).  The 
coding  of  the  present  EVAL  is  entirely  new  and  represents  a 
different  approach  to  many  (although  not  all)  of  the 
computations . 

4 • 2 USER  OPTIONS  AND  FLEXIBILITY 

The  EVAL  program  is  quite  flexible;,  the  user  has 
the  necessary  control  of  program  flow  through  input  options . 
The  options  are  set  through  punch  card  inputs  read  at  the 
start  of  processing.  The  various  inputs  (along  with  their 
formats)  are  listed  in  Table  4-1  and  described  below. 

First,  the  user  determines  which  target  types  are 
to  be  considered  during  the  EVAL  run.  Types  are  numbered 
from  1 to  30;  the  actual  grouping  of  targets  into  types  is, 
of  course,  determined  at  the  time  of  SCENARIO  execution,  but 


* See  Appendix  F,  AIRS,  Vol.  I,  for  definition  of  the  beta 
measure . 
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the  user  has  the  option  of  excluding  any  particular  target 
type  from  the  processing  done  by  EVAL.  This  option  is 
exercised  by  setting  the  input  variable  IHOLD(J),  J = 1,30 
The  input  values  are  simply  IIIOLD(J)  = 1 if  type  number  J 
should  be  included  in  the  processing,  and  IHOLD(J)  = 0 if 
it  is  to  be  excluded. 


Second,  the  user  specifies  a maximum  offset 
distance . If  the  distance  from  the  aircraft  to  a particular 
target  sighting  is  greater  than  this  maximum,  the  sighting 
is  not  considered  by  EVAL.  (This  action  does  not  affect 
other  sightings  of  the  same  target,  or  other  targets  of  the 
same  type.)  The  input  variable  for  this  option  is  AMAXOF. 

Last,  the  user  may  define  up  to  five  target  groups , 
where  a group,  as  explained  above,  is  a collection  of  target 
types.*  The  input  variables  are: 


where 

and 


NUM(I),  1 = 1,5 

NTY (I , N)  I = 1,  5;  N = 1 , NUM(I) 

NUM(I)  is  the  number  of  target  types  in  group  I, 


NTY  (I,  N)  is  the  type  number  (1  to  30)  of  the  Nth 
type  of  the  Ith  group. 

Summary  measures  are  computed  and  printed  for  each  of  the 
five  groups. 


* Unlike  the  previous  EVAL  routine,  groups  need  not  be 
disjoint.  That  is,  a particular  type  may  belong  to  more 
than  one  group. 


1 

t 

i 


\ 

i 


i 

% 

s 

s 

« 

i 

■f 

<* 


K 

*. 


\ 

r% 

L 


V 

V 


The  flexibility  of  EVAL  can  now  be  illustrated.  As 
mentioned  before,  the  EVAL  routine  is  coded  as  a loop 
around  the  processing  phase,  so  the  user  may  evaluate  a 
particular  emulation  result  (as  recorded  on  the  Target/ 
Sensor  Output  Tape)  in  many  ways  by  stacking,  one  behind 
the  other,  sets  of  user  option  parameter  cards.  By  using 
AMAXOF,  the  user  can  determine  the  effects  of  varying  maxi- 
mum offset  distances  on  system  performance  in  terms  of 
target  types  and  groups.  Moreover,  the  user  can  evaluate 
system  performance  with  respect  to  any  collection  of  target 
types.  For  a quick  look  at  system  performance  (i.e.,  total 
system  detectability,  identi f iabi li ty , and  locali zabi li ty 
at  each  processing  level)  , the  user  groups  the  relevant 
types  together  through  the  use  of  NTY  and  then  reads  the 
target  group  output  tables.  For  a more  detailed  break  out 
he  simply  "turns  off"  the  target  types  he  is  not  interested 
in  (through  the  use  of  the  IHOLD  array)  and  receives  the 
total  EVAL  output  with  only  the  right  types  considered.* 

4 • 3 INPUT 

EVAL  uses  both  tape  and  punch  card  inputs.  The 
tape  input  consists  of  the  Target/Sensor  Output  Tape, 
previously  referenced.  The  tape  consists  of  records  con- 
taining the  necessary  data,  one  record  per  target  sighting. 
EVAL  requires  that  the  Target/Sensor  tape  be  mounted  on 
logical  unit  2.  The  punch  card  input  consists  of  8 cards; 
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* The  user  option  inputs  NONPRE  and  IPREPL  used  in  the 
previous  EVAL  to  achieve  preplanned  vs.  non-prcplanncd 
flexibility  have  been  eliminated. 
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as  many  packs  of  8 data  cards  as  the  number  of  desired 
loops  through  the  program  are  needed'.  A blank  card  should 
be  placed  behind  the  last  pack  of  8 to  end  processing; 
the  CDC  operating  system  requires  that  an  end-of-file  card 
follow  this  blank  one.  Seven  of  the  8 input  cards  are 
concerned  with  the  user  option  variables  described  above; 
the  other  is  a card  containing  an  alphanumeric  header.  The 
header  is  printed  out  in  addition  to  the  header  obtained 
from  the  EXECUTIVE  tape;  thus,  the  user  should  use  a differ- 
ent header  card  in  each  pack  of  8 ( if  he  uses  more  than  one 

pack  at  a time)  so  that  the  EXECUTIVE  header  serves  as  an 

overall  job  title,  and  the  header  cards  serve  as  run  sub- 
titles. The  exact  form  of  each  input  card  as  well  as  their 

order  is  shown  in  Table  4-1. 

4.4  PROCESSING 

The  EVAL  routine  produces  output  via  the  following 
series  of  steps  : 

Step  0 : 

The  user  input  options  are  read  from  punch  cards. 

Step  1 : 

The  Target/Sensor  Output  Tape  is  positioned,  and 

all  variables  and  accumulating  arrays  are  initial  - 

i zed . 

Step  2: 

The  first  physical  record  of  a target  sighting 

record  is  read  from  the  tape.* 


* A target  sighting  record  on  the  Target/Sensor  Output  Tape 
which  contains  all  the  needed  data  on  the  sighting,  consists 
of  20  physical  tape  records;'  see  reference  to  the  tape 

given  earlier. 


4-6 


m * * m "m  9J>  ' J'  ’/  ■** , **.  * ■ ^ - 


*■_ . 


Step  3: 


Next 


A check  of  IHOLD  and  AM AX OF  is  made  against  the 
appropriate  target  sighting  data.  If  the  target 
sighting  is  not  to  be  considered,  all  further 
physical  records  of  the  sighting  are  skipped,  the 
tape  is  positioned  at  the  beginning  of  the  next 
sighting,  and  control  returns  to  Step  2.  System 
performance  data  is  not  affected.  If,  on  the  other 
hand,  the  sighting  is  to  be  considered,  processing 
continues  with  Step  4. 

Step  4 : 

The  target  sighting  data  broken  out  by  sensors 
is  read  from  the  tape.*  A check  is  made  to  deter- 
mine whether  at  least  one  sensor  had  the  target 
sighting  within  its  field  of  view  (FOV) . If  not, 
all  further  records  are  skipped  (as  in  Step  3)  and 
control  returns  to  Step  2 for  the  processing  of  a 
new  sighting;  if  so,  the  data  is  assimilated  into 
the  proper  accumulating  arrays.  This  consists 
of  accumulating  into  arrays  indexed  by  the  sighting's 
target  type  number  (which  will  produce  the  effec- 
tiveness measures  of  types  in  sighting  units)  and 
by  the  sighting's  target  number  (which  will  allow 
the  computation  of  target  measures  conjoined  from 
individual  sightings,  and  eventually  produce  type 
measures  in  target  units)  . 

Step  5: 

Overall  system  performance  data  which  has  been 


4 physical  records. 
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integrated  over  sensors  by  EXECUTIVE  are  next 
read  in  by  EVAL  and  accumulated  into  arrays  as  in 
Step  4.*  Control  then  returns  to  Step  2 to  begin 
the  processing  of  the  next  target  sighting.  If, 
however,  no  more  sightings  remain  on  the  tape, 
control  goes  to  Step  6. 


With  the  processing  of  target  sightings  complete, 
the  accumulating  arrays  indexed  by  target  number 
are  processed  to  form  target  measures  - that  is, 


>V- 
> .** 


measures  conjoined  over  sightings  of  a particular 
target . 


Step  7 : 

The  target  measures  are  accumulated  into  types. 

At  this  point  then  there  are  two  kinds  of  arrays 
indexed  on  type  --  those  containing  target  unit 
data  and  those  with  sighting  unit  data.  The  latter 
where  formed  in  Step  4,  the  former  in  this  step. 

Step  8: 

All  arrays  of  both  kinds  indexed  on  type  are  processed 
to  form  final  type  measures  in  both  target  and 
sighting  units. 

Sten  9 : 

The  type  measures  are  then  integrated  to  form  over- 
all system  measures  in  both  units.  Note  that  although 


type  measures  are  being  integrated,  the  integration 
is  done  in  such  a way  that  the  system  measures  are 
in  fact  taken  over  targets  and  sightings,  not  over 
targets  and  sightings  first  averaged  into  types. 

For  ease  of  output  and  internal  indexing,  these 
final  measures  are  given  the  label  "Target  Type  31" 
(there  are  only  30  real  types). 

Step  1 0 : 

In  accordance  with  the  user  defined  NTY  array, 
the  target  and  sighting  unit  measures  of  the 
appropriate  types  are  integrated  to  form  the 
target  group  measures.  This  integration  is  per- 
formed similarly  to  that  of  Step  9;  this  results 
in  properly  averaged  measures. 

Step  11: 

The  computed  measures  are  outputted. 

Step  12: 

If  the  next  card  in  the  user  punch  card  data  deck 
is  blank,  the  EVAL  program  ends  and  control  is 
passed  to  the  CDC  monitor.  Otherwise  control  goes 
to  Step  0 for  another  processing  of  the  same  tape 
with  different  user  options. 

The  processing  described  in  Steps  0 - 12  is  flow- 
charted in  Figure  4-1.  During  this  processing,  the  EVAL 
routine  uses  the  disk  for  intermediate  storage  and  employs 
the  NAVAIRDEVCEN  disk  software.  All  disk  communications  are 
handled  by  two  EVAL  subroutines  entitled  GET  and  l’UT.  The 
former  is  a two  argument  function  subprogram,  and  the  latter 
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a subroutine  subprogram.  GET  and  PUT  perform  disk  allocation, 
maintain  a disk  memory  map  through  the  use  of  COMMON  variables, 
compute  disk  optimization,  and  make  all  calls  to  disk  soft- 
ware. EVAL  requires  the  disk  to  be  defined  as  logical  unit  14. 
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4.5  COMPUTATIONS 

The  following  quantities  are  computed  by  the  above 
described  processing.  The  computations,  except  where  indicated, 
are  straight-forward  averaging  or  summing  of  data  contained 
on  the  Target/Sensor  Output-  Taoe. 

Data  broken  out  by  sensor  and  type: 

(la)  The  number  of  times  targets  of  the  particular 
type  fell  within  the  FOV  of  the  sensor. 

(lb)  The  number  of  times  target  sightings  of  the 
particular  type  fell  within  the  FOV  of  the 
sensor . 

(2a)  The  mean  probability  that  the  targets  were 
not  masked  by  terrain. 

(2b)  The  mean  probability  that  the  target  sightings 
were  not  masked  by  terrain. 

(5a)  The  mean  probability  that  the  particular 

sensor  was  not  non-operational  when  viewing 
the  targets. 


The  mean  probability  that  the  particular 
sensor  was  not  non-operational  for  the 
target  sightings. 


The  mean  probability  that  the  targets  were 
not  obscured  by  clouds. 

The  mean  probability  that  the  target 
sightings  were  not  obsured  by  clouds. 

The  expected  number  of  times  the  targets 
were  detected  at  levels  1,  2,  3,  and  4. 

The  expected  number  of  times  target  sighting 
detections,  were  made  at  levels  1,  2,  3,  and  4 

The  expected  number  of  times  the  targets  were 
identified  at  levels  1,  2,  3,  and  4. 

The  expected  number  of  times  target  sighting 
identifications  were  made  at  levels  1,  2, 

3,  and  4. 

The  mean  conditional  CEP  in  feet  of  the 
targets . 

The  mean  conditional  CEP  in  feet  of  the 
target  si ghtings . 

The  beta  measure  of  sensor  redundancy  for 
targets . 

The  beta  measure  of  sensor  redundancy  for 
target  sightings. 


V ILVV5>>.\Vf 


Concerning  items  (2a),  (3a),  and  (4a),  the  probability 
for  a single  target  (i.e.  one  of  the  valves  averaged  to  form 
the  indicated  mean)  is  defined  as  the  probability  of  success 
on  all  sightings  of  the  target.  Note  that  items  (2a),  (2b), 

(3a)  , (3b)  , (4a)  , and  (4b)  are  means  of  the  probabilities  of 
non-independent  events;  care  should  be  taken  in  using  them. 
Concerning  items  (5a)  and  (6a),  an  expected  detection  or 
identification  of  a single  target  is  defined  as  the  expected 
detection  or  identification  on  at  least  one  sighting  of  the 
target.  Concerning  items  (.7a)  and  (7b),  the  CEP  combining 
algorithm  (Paragraph  3.3.4,  AIRS,  Vol.  I)  is  used  on  sighting 
CEP's  to  produce  a target  CEP.*  Information  on  items  (8a) 
and  (8b)  is  found  in  Appendix  F,  AIRS,  Vol.  I;  it  is  computed 
similarly  for  targets  and  sightings.  Since  the  Target/Sensor 
Output  Tape  contains  only  individual  sighting  records,  the 
conjoining  over  sightings  into  single  target  measures  re- 
quired by  (2a),  (3a),  (4a),  (5a),  and  (7a)  is  done  by  EVAL 
(Processing  Step  6) . 

Data  broken  out  by  type: 

(9)  Total  (i.e.,  conjoined  over  the  sensors) 
target  unit  and  sighting  unit  mean 


* An  investigation  was  made  to  determine  how  the  approximation 
embodied  in  the  CEP  combining  algorithm  fared  when  applied  to 
the  small  number  of  sightings  of  a particular  target,  and  also 
the  feasability  of  using  exact  solutions  in  this  particular 
case.  It  was  concluded  that  although  exact  solutions  could 
be  employed,  it  would  not  be  worth  the  effort  because  the 
approximation  appeared  quite  sound.  The  path ol ogi cal  cases 
turn  out  to  be  not  those  with  a small  number  of  measurements, 
but  rather  those  with  grossly  unequal  measurements. 
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detectability  at  each  of  the  four  levels. 


(10)  Total  target  unit  and  sighting  unit  mean 
identif ia’oility  at  each  of  the  four  levels. 

(11)  Relative  target  and  sighting  mean  CEP  at 
each  of  the  four  levels. 

(12)  Total  target  and  sighting  mean  CEP  at  each 
of  the  four  levels. 

Individual  sighting  total  measures  are  taken  from 
the  Tape  (Processing  Step  5) ; individual  target  total  measures 
are  formed  in  the  same  way  as  (5a),  (6a),  and  (7a). 

4.6  OUTPUT 

The  printed  reports  generated  by  a single  run  of 
EVAL  are,  in  order: 

(a)  A header  page  containing  job  title,  run  title, 
the  numbers  of  the  target  types  being  considered, 
and  the  maximum  offset  distance. 

(b)  Pages  of  tables  containing  the  target  type 
statist! cs  for  each  type  considered  (one  page  per 
type).  The  statistics  consist  of  data  items  (la) 
through  (12)*  and  the  additional  information: 

. The  number  of  targets  being  considered 
in  the  target  type. 

. The  number  of  target  sightings  being 
considered  in  the  target  type. 


* At  the  present  time,  item  (8)  is  not  printed. 
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. The  number  of  targets  of  the  type  falling 
into  the  FOV  of  at  least  one  sensor. 

The  number  of  target  sightings  of  the 
type  falling  into  the  FOV  of  at  least 
one  sensor. 


(c)  A one  page  table  containing  the  overall  sys tern 
performance  measures  (target  type  #31)  in  the  same 
format  as  the  individual  type  tables,  with  the 
exception  of  data  item  (8) . 


(d)  Two  pages  containing  the  overall  system  re  - 
dundancy  measures  in  both  sighting  and  target  units 
(data  item  (8) ) . 


(e)  Two  pages  containing 
which  include  items  (10), 
These  pages  also  list  the 
each  target  group. 


the  group  statis  tics , 
(11) , and  (12)  only, 
target  types  included  in 


4.7  MISCELLANEOUS 


The  EVAL  routine  is  coded  in  FORTRAN,  wel 1 - commented , 
and  quite  readable.  To  aid  in  further  understanding  the  code, 
the  implementers  have  used  certain  prefixes  and  suffixes  to 
designate  certain  kinds  of  variables: 


Code 


Meaning 


S 


T 


A variable  being  used  to 
accumulate  sighting  unit 
dat  a . 

A variable  being  used  to 
conjoin  sighting  data 
into  target  measures. 


Code  Me  aji  1 n g 

F A variable  being  used  to 

integrate  over  targets  to 
form  target  unit  data. 

A A variable  containing  data 

integrated  over  sensors. 

G A variable  containing  data 

integrated  into  groups. 


Moreover,  when  used  as  array  subscripts,  certain  integer 
constants  have  reserved  meanings:  IT  always  refers  to  target 

type,  L to  level  number,  N to  group  number,  I to  target 
number,  and  J to  sensor  number. 

Furthermore,  the  variables  have  been  made  as  meumonic 
as  possible.  For  example  then,  the  variable  APDS  (IT,L)  would 
be  the  total  (integrated  over  sensors)  probability  of  de- 
tection in  sighting  units  for  level  L and  target  type  IT. 
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SECTION  V 


SENSOR  MODEL  IMPROVEMENTS 


5.1 


INTRODUCTION 


Analytics  was  requested  to  make  improvements  to 
the  various  sensor  models  reported  in  the  AIRS,  Vol.  I. 
These  improvements  are  discussed  in  this  section. 


•S*//VJ 


Briefly,  the  following  improvements  have  been  made 
to  the  sensor  models: 


IR  and  FL1R:  These  models  have  been  improved 

by  allowing  for  image  motion  (i.e., 
improper  V/H  compensation)  with  re- 
spect to  sensor  mot  ion  to  degrade 


detectabilities  and  ident i f i abi litics . 
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SLR:  This  model  has  been  improved  by  allowing 

for  the  use  of  either  synthetic  or  real 
aperture  processing.  Also  considered  is 
the  fact  that  only  a certain  range 
interval  is  recorded  on  film  (i.e.,  only 
a certain  part  of  the  actual  field  of  view) . 

MTISLR  and  MTIFLR:  These  sensor  models  have 

been  improved  by  taking  into  account 
additional  equipment  parameters  such  as 
(1)  the  actual  sensor  response  as  a func- 
tion of  the  doppler  shift  of  signals 
received  from  a target,  (2)  the  number  of 
delay  lines,  and  (3)  other  such  parameters. 
Furthermore,  other  effects  of  the  environ- 
ment such  as  the  direction  and  speed  of 
rainfall  are  considered. 

ECM:  This  model  now  explicitly  considers  the 

possibility  that  signals  will  not  be 
detected  because  of  the  density  of  the 
EM  environment. 

PHOTO:  The  following  improvements  have  been 

made  and  incorporated  into  this  model 

image  motion  may  now  degrade  camera 
performance,  especially  the  identi- 
fiability  of  targets. 

the  effects  of  atmospheric  transmis- 
sivity and  backscatter  on  night 
photography  are  now  considered. 

for  night  photography,  the  illuminator 
source,  power,  directivity,  distance 
from  target  and  camera,  and  other 
associated  parameters  are  now  considered. 

a more  realistic  computation  for  contrast 
during  daylight  hours. 

independence  of  cloud  cover  probability 
on  successive  frames  lias  now  been  discard- 
ed; dependence  is  considered. 


These  improvements  substantially  improve  the  models  and  extend 
the  capabilities  of  the  AIRS  Model  by  allowing  cons iderat ion 
of  parameters  heretofore  either  not  treated  at  all  or  treated 
only  superficially.  The  IR  and  FLIR  improvements  are  dis- 
cussed in  Paragraph  5.2  of  this  report,  and  SLR  improvements 
are  discussed  in  Paragraph  5.3.  The  MTISLR  and  MTIFLR 
improvements  are  the  topics  of  Paragraphs  5.4  and  5.5,  re- 
spectively. Improvements  in  the  ECM  model  arc  discussed  in 
Paragraph  5.6  and  those  in  the  PHOTO  model  in  Paragraph  5.7. 

5.2  IR  AND  FUR 

5.2.1  Introduction 

An  IR  or  FLIR  sensor  will  occasionally  miss  de- 
tecting a target  because  the  sensor  has  not  been  properly 
adjusted  to  take  into  account  the  interrelated  effects  of 
the  aircraft  velocity  and  height  above  ground. 

Each  time  the  sensor's  mirrors  rotate,  a path  is 
traced  out  on  the  ground  (the  IR  can,  at  least  potentially, 
"see"  a target  that  falls  within  this  path) . The  thickness 
of  this  path  on  the  ground  is  a function  of  three  parameters: 
the  aircraft  height  above  ground,  the  ground  offset  distance 
from  the  aircraft  vertical,  and  the  angular  beam  spread  of 
the  detector.  A short  time  later,  as  the  mirrors  again 
rotate,  a second  path  is  traced  on  the  ground.  Since  the 
aircraft  lias  moved  forward  during  the  time  between  scans, 
the  second  scan  will  not  sweep  over  the  same  ground  as  the 
first.  If  the  aircraft  velocity  is  too  great,  the  two  paths 
do  not  overlap  and  there  is  a gap  between  successive  scans 
which  the  IR  does  not  "see".  Figure  5-1  illustrates  nic- 
torially  the  effect  of  the  Y/ii  ratio. 
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The  following  paragraph  will  explore  how  this 
velocity/height/time-bctween-scans  mismatch  affects  the  per- 
formances of  the  IR  and  FLIR  sensor.  Specif ically , Paragraph 
5.2.2  will  define  terms  to  be  used  and  develop  geometric 
equations  which  are  applicable  to  both  the  IR  and  FLIR 
sensors.  Paragraph  5.2.3  will  then  develop  for  the  IR, 
equations  which  yield  the  probability  that  a given  target 
does  not  fall  into  a scan  gap  and  shows  how  this  probability 
affects  the  detectability  of  that  particular  target. 

Paragraph  5.2.4  presents  an  analogous  treatment  for  the  FLIR. 


5.2.2 


Definition  of  Terms.  Basic  Geometric  Equations 


First,  an  expression  will  be  developed  for  the 
aircraft  location  at  the  time  a point  P on  the  ground  is 
scanned  by  the  IR  or  FLIR  sensors.  The  most  general  case, 
slewed  FLIR,  is  considered. 

First  consider  the 
axis  of  the  sensor  to  be  slewed 
through  an  angle  \p  in  the 
horizontal  plane,  \p  being  positive 
to  the  right.  In  Figure  5-2,  the  A Xg 

subscripts  A denote  aircraft 

axes  and  the  subscripts  B denote  / 

the  coordinate  axes  after  rotation 
through  . The  aircraft  is  at  1 

the  common  origin  q.  The  common  / 

Z^  or  Zg  axis  is  positive  down-  / 

ward.  X! 


Ai  rtfraf  r 


Figure  5-2 
Plane  View 


Figure  5-3  shows  aircraft 

a vertical  plane  through  \ Xj-. 

the  Xg  axis.  In  ordi-  s' 

nary  IR,  the  mirror  axis  N.  s'  XO 

of  rotation  is  horizon-  ) 

tal,  about  Xg.  In  FLIR,  s'  ^ ^ 

the  axis  is  tilted  up  \ \ 

through  an  angle  8 . The  \ ^ 

subscript  C denotes  the  ^ 

tilted  axis.  The  Yc  axis  g \ / 

D \ V (X 

is  common  to  the  coordi-  * 

nate  systems  before  and  - . ...  P \ 

after  tilting.  It  is  . ' . ' .• 

assumed  that  the  mirrors  vl  - Z0 

ZA,ZB  C Figure 

look  out  at  right-angles  Vertic 

from  the  axis  of  rotation,  rather  than  at  some  arbitrary  ^ 
angle  a as  shown  dotted  in  Figure  5-3. 

The  mirrors  rotate  through  an  angle  <J>  = tot  about 
the  X^  axis.  As  they  do,  the  intersection,  between  OZ^, 
and  the  ground  plane  traces  a path  in  the  ground  plane.  This 
is  the  path  of  the  point  in  the  center  of  the  "beam”.  P 

o 

in  Figure  5-3  is  part  of  the  path.  Let  P denote  an  arbi- 
trary point  on  this  path. 


Letting : 


where 


H = H, 


^AC  = heiSht  of  aircraft  above  ground 
- height  of  target  above  ground 
resolutions  through  the  angles  , , and  give  for  the 
point  P coordinates: 

X = H(sin  0 cos  <J>  cos  ip  - sin  <p  sin  ip)  / cos  0 cos  <J> 

Y = H(sin  G cos  <j>  sin  vJj  + sin  <f>  cos  ip~)  / cos  G cos  4> 

Z = 11, 

where  X,Y,  and  Z are  in  the  aircraft  coordinate  system, 

X being  along  the  flight  line. 
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Ordinarily,  the  target  position  with  respect  to 
the  flight  line,  i.e.,  its  Y-coordinate , is  known,  and  it  is 
desired  to  find  the  X-coordinate  at  the  time  of  sighting. 

The  value  of  X is  found  by  eliminating  <J>  between  equations 
(5.1a)  and  (5.1b);  the  result  for  a point  target  at  the  center 
of  the  "beam"  is: 

X = (H  tan  9 - Y sin  ^)  / cos  (5.2) 

The  range  at  the  time  the  target  is  scanned  is: 

R(Y)  - (H2  + X2  + Y2)^  (5.3) 

where  X is  given  by  equation  (5.2). 
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5.2.3 


Ordinary  IR 


In  ordinary  IR,  equations  (5.1a),  (5.1b)  and  (5.1c) 
are  simplified  because  tp  and  8 are  both  zero,  yielding: 

X = 0 ; Y = H tan  <j> 


R(Y)  = (H2  + Y2)' 


(5.4) 


There  are  two  common  scanning  modes.  In  one  the  detectors 
scan  in  unison:  <f>  has  the  same  value  for  all  detectors  at 

any  instant.  In  the  other  mode  the  detectors  are  evenly 
spaced  around  the  cycle.  If  the  angular  rotation  uj  does 
not  match  the  aircraft  forward  velocity,  the  overlap  (if  u 
too  large)  or  the  gap  (if  w too  small)  is  evenly  distributed 
between  adjacent  detectors  if  they  are  staggered,  and  is 
lumped  if  the  detectors  scan  in  unison.  The  staggered  method 
is  chosen  as  more  nearly  representative  of  today's  IR.  There- 
fore,  the  j detector  "sees"  at  the  angle  <t>  ^ , given  by: 

4>.  = ut  + 2tt j / 

where  is  the  number,  w the  angular  frequency,  and  t 

the  time  between  scans  of  the  detector.  Since  the  staggered 
scan  mode  is  used,  it  is  unnecessary  to  assign  different 
angles  of  a (sec  Figure  5-3)  to  the  different 
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detectors;  a is  zero  for  all  of  them.  If  V is  the  air- 
craft forward  velocity,  the  beam  centers  for  successive 
detector-scans  are  separated  in  the  X-direction  by  the 
distance  in  the  ground  plane, 

AX  = VT/Xd  = 2ttV/uNd  (5.5) 

where  T - 2tt/u  is  the  period  of  rotation  of  the  mirrors. 

The  beam  width  in  the  ground  plane,  for  a point 
target  at  P,  is 

SR(Y)  = 6 (H2  + \2)h 

where  6 is  the  beam  width  in  radians.  For  point  targets 
at  distance  Y off  the  flight  path  the  difference  U,  where 

U = 6R(Y)  - AX  (5.6) 

is  the  amount  of  overlap  between  successive  detector  scans 
if  U is  positive;  U is  the  gap  between  successive  detector 
scans  if  it  is  negative.  In  the  event  that  a gap  occurs,  the 
probability  that  a point  target  with  coordinate  Y is  in 
the  gap  between  scans  is 

(AX  - 5R (Y) ) / AX  (5.7) 

and  the  probability  that  the  target  is  scanned  and  not  missed, 


PN  * 1S 


PN  = Min 


{ 1,  6R(Y)  / AX  } 


(5.8) 


Pm  is  next  adjusted  to  account  for  targets  not 


N 

being  "point  targets": 


£ 


% 


PN,  = Min 


{ 1,  ( R (Y)  + \ ) / AX  J (5.9) 
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where  L is  the  length  of  the  target  in  the  X direction 
and  Y is  the  Y-coordinate  for  the  center  of  the  target. 

This  probability  is  used  multiplicatively  in  the 
equations  for  total  detectability.  Specifically,  for  pro- 
cessing level  1,  the  total  detectability,  PD,  of  a given 
target  becomes: 

(5.10) 


PD  = PDT  (Px)  (P2)  (Kl)  (Pn) 


where : 


PDT 

P, 


the  conditional  detectability  for  level  1 
the  terrain  non-shadowing  probability 
P2  = the  equipment  up  probability 

= the  decision  (0  = no,  1 = yes)  of  whether 
the  Ik  sensor  operates  in  level  1. 

Similarly,  the  total  detectability  for  levels  2,3,  and  4 
of  processing  becomes 


PD 


= PCD  (Px)  (P2)  (Kx)  (PN) 


(5.11) 


where : 


PCD  = the  conditional  detectability  for 
levels  2,  5,  and  4 

= the  decision  (0  = no,  1 = yes)  of  whether 
the  IR  sensor  operates  in  the  level  (2,3, 
or  4)  being  considered  and  all  other  terms 
are  as  defined  above. 

FLIP  --  Forward  Looking  Infrared 

There  are  two  modes  in  which  FLIP  is  operated:  a 
not-slewed  mode,  and  a slewed  mode.  For  the  not-slewed 
mode,  the  sensor  is  always  directed  forward  and  downward. 
In  the  slewed  mode,  the  sensor  is  pointed  toward  a known 
target,  typically  off  the  flight  line.  Those  two  modes 


5.2.4 
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must  be  modeled  separately.  Although  the  computation  of  the 
probabilities  of  detection  for  each  mode  are  somewhat  modi- 
fied, the  method  of  combining  the  modes  is  unchanged  from  that 
reported  in  Paragraph  3. 2. 8.1  of  AIRS,  Vol.  I. 

5. 2. 4.1  FLIR,  Not-slewed.  For  FLIR  not-slewed,  is  0; 
equations  (5.1a),  (5.1b)  and  (5.1c)  become: 

X = H tan  0 
Y = H tan  <j>  / cos  0 
1 = H 

and  equation  (5.3)  becomes 

R(Y)  = (Y2  + H2(l  + tan2  0 ) ) ^ (5.12) 

The  general  reasoning  that  led  to  equation  (5.7)  may  be 
applied  here  also  but  instead  of  the  beam  width,  the  ground- 
projected  beam  width  must  be  used: 

R(Y)  / sin  0 (5.13) 

And  in  place  of  L/2,  half  of  the  ground-projected  length 
must  be  used 


L/2  for  flat  targets 
(L/2)  / sin  0 for  spherical  targets, 
so  that  the  probability,  PjN[,  that  the  target  does  not  fall 
in  a scan  gap  is: 


PN  = 


Min  {1,  (SR (Y)  + L sin  0/2)  / AX  sin  0}  TP*oet< 
Min  {1,  (SR (Y)  + L/2)  / AX  sin  0}  Argots1 


(5.1 
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Although  equations  are  given  for  both  flat  and 
spherical  targets,  it  was  decided  to  model  all  targets  as 
spherical.  Equations  (5.10)  and  (5.11)  are  then  used  to 
compute  the  probability  the  FLIR  detected  the  target  in  the 
not-slewed  mode. 


5. 2. 4. 2 FLIR  Slewed.  For  slewed  FLIR,  equations  (5.1)  and 
(5.3)  are  used  as  written.  The  point  P still  traces  a 
straight  line,  but  the  slope  of  the  line  is  no  longer  zero. 
It  is  given  by 


dX  . . ' 

37  = - tan  ip 

Therefore,  the  distance  be- 
tween beam  centers  is  not 
AX  as  given  by  equation 
(5.5),  but  is  actually: 


X flight  path 

trace  for 

! next  detector 


trace  of  P 


AX  cos  i p 


AX  j”  cos  \p  | 

and  this  is  the  distance 
that  must  be  used  in  place 
of  AX  in  equation  (5.13), 


H tan0/cosi|) 


Then 


Figure  5-4 
Ground  Plane  View 


(<5R  (Y)  + L sin  (f ) ) 


| ( jeos^  | ) AX  sin  0 ) 

PN  = i 

(6R(Y)  + (1-)) 

Min  1 , £ 

( |cosi/|)AX  sin  0 

for  flat  and  for  spherical  targets,  respectively. 


(5.15) 


(5.16) 
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As  in  the  ease  of  FLIR  not-slcwed,  a spherical 
target  is  assumed.  Equations  (5.10)  and  (5.11)  can  now  be 
used  to  give  the  probability  of  detection  for  FLIR  in  (1) 
level  1,  and  (2)  levels  2,3,  and  4,  respectively. 


5.3  SLR 

The  SLR  model  has  been  extended  to  allow  either 
synthetic  or  real  aperture  radars  to  be  modeled.  Heretofore, 
only  synthetic  aperture  radar  was  modeled,  and  cross-range 
resolution  was  held  to  be  the  same  in  all  four  processing 
levels.  Now,  however,  if  a synthetic  aperture  radar  is 
being  modeled,  consideration  is  given  to  the  fact  that  the 
cross-range  resolution  attainable  in  level  1 is  different 
from  that  attainable  in  the  other  processing  levels. 


5.3.1  Synthetic  Aperture 

If  the  user  specifies  synthetic  aperture,  he  must 
also  stipulate  (1)  the  cross-range  resolution  attainable  in 
level  1,  and  (2)  the  cross-range  resolution  attainable  in 
levels  2,  3,  and  4.  The  attainable  resolution  in  level  1 is, 
of  course,  less  than  that  of  the  other  levels;  this  is  due 
primarily  to  lack  of  time  to  adequately  process  the  signal 
returns . 


Equation  (3.36)  of  AIRS 
ground  resolution  cell  size,  A 


Vol.  I,  which  gives 
may  be  rewritten  as: 


the 


Ar  " ^1  ( ^p2  "+  h (2R+h)  - p ) 
Ar  ■^2^'/p2+  h(2R+h)  - p ^ 


Level  1. 
Level  2,3  or  4. 


1 1 


where 


1WI+J  • v 


5.3.2 


= the  cross-range  resolution  attainable  in  level  1 

= the  cross-range  resolution  attainable  in  levels 
2 , 3 or  4 . 

= ground  range  to  target 
= range  resolution 
= slant  range  to  target 


Real  Aperture 


If  real  aperture  (sometimes  called  brute  force) 
radar  is  being  modeled,  the  resolution  attainable  in  all 
four  processing  levels  is  a direct  function  of  both  the  beam 
width  and  slant  range  to  the  target.  In  particular,  the  cross 
range  resolution  for  real  aperture  is  given  by  $R;  where 
3 is  the  half-power  beam  width  and  R is  the  slant  range  to 
the  target.  The  half-power  beam  width  is  selected  because  it 
is  consistent  with  what  is  usually  meant  by  the  term  "beam 
width".  Thus,  rewriting  equation  (3.36)  of  AIRS,  Vol.  I for 
the  real  aperture  case: 


Ar  - BR  ( 


p2  + h(2R+h)  - p 


The  symbols  are  defined  in  Paragraphs  5.3.1  and  5.3.2 


5.3.3  Detectability,  Idcntif iabil i ty , and  Local inability 

A check  was  added  in  the  SLR  model  to  make  certain 
that  the  target  fell  within  the  maximum  and  minimum  slant 
ranges  recorded  on  film.  If  the  target  did  not  fall  within 
these  limits,  detectability  and  identifiability  for  levels 
2,  3,  and  4 were  set  to  zero,  as  these  levels  require  the 
targets'  appearance  on  film.  Also,  if  real  aperature  is  used 
A (cqn . 3.74  of  AIRS  Vol.  I)  is  set  to  zero.  With  these  ex- 
ceptions, all  equations  given  in  AIRS,  I bl.  I for  computing 
detectability , ident ifiabi 1 ity,  and  the  CLP  remain  unchanged  - A 


A 
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The  detection  of  moving  targets  by  side  looking 
radar  (MTISLR  model)  has  been  refined  by  the  more  precise 
calculation  of  previously  approximated  equipment  and  environ- 
mental variables.  The  calculations  permit  the  model  to  take 
into  account  the  effects  of  radar  instabilities,  aircraft 
speed,  precipitation  particle  velocity,  and  wind  clutter. 

These  refinements  are  available  at  the  cost  of  the  user  having 
to  provide  about  twice  as  many  input  parameters  as  were  needed 
for  the  simpler  model. 

5.4.1  Outline  of  Revised  Computations 

With  the  aid  of  two  assumptions  (one  dealing  with 
the  elimination  of  the  effect  of  aircraft  motion  from  the  Dop- 
pler frequency  shifts  and  the  other  dealing  with  the  nature 
of  the  MTI  filter),  which  are  developed  below,  the  output  of 
the  radar  processing  can  be  computed  for  the  return  from  any 
resolution  unit.  For  any  particular  resolution  cell,  contain- 
ing a target,  the  output  power  is  called  S.  If  the  target  is 
replaced  by  background  and  the  output  power  computed,  the  re- 
sult is  called  C.  Also  the  response  of  the  filter  to  the 
radar  noise  can  be  computed,  if  the  input  noise  character- 
istics are  known:  call  the  output  noise  N.  Then  the  ratio 
S/(C+N) , can  be  used  as  the  argument  for  the  cubic  polynomial 
given  in  the  AIRS  report,  Vol.  I.,  as  equation  (3.103),  p.5-50, 
which  gives  the  probability  of  detection  P^. 

5.4.1. 1 Doppler  Shift  of  the  Echo.  For  any  target,  or  for 
any  area-element  of  precipitation  if  any,  equations  (5.17), 
(5.18),  and  (5.19)  provide  the  means  for  computing  the 
Doppler  shift  of  the  echo. 
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Let  (XA,  Y^,  Z^)  be  the  position  of  the  aircraft 
and  = dX^/dt  its  velocity,  assumed  in  the  X-direction. 
Let  (X, Y, Z)  be  the  coordinates  and  (V^,  Vy , Vz)  the 
velocity- vector  components  of  an  object  from  which  a radar 
echo  is  received.  Then  the  radar  range  is 

R = [(X-XA)2  + (Y-Ya)2  + cz-zA)  2J%  (5.17; 

and  the  rate  of  change  with  time  is: 

at  “ * = <VV  * <y-yA>  VV  + (Z-XA>  Vzl 


(5. IS) 

The  Doppler  frequency  shift  in  the  received  echo  is: 

fD  = -2C§)  fR  (5.19) 

Where  c is  the  speed  of  light  and  fR  is  the  radar 
transmitted  frequency. 

Since  the  power  in  the  echo  is  expressible  for  each 
such  element  via  the  computation  of  radar  reflectivities, 
the  power  spectrum  of  the  echo  can  be  computed  for  any  given 
resolution  cell  by  summing  over  all  the  reflecting  elements 
in  the  cell. 

After  the  spectrum  of  the  echo  has  been  computed, 
the  next  step  is  to  compute  the  output  of  the  radar  filter 
in  response  to  an  input  with  that  spectrum.  In  order  to 
compute  the  filter  output  two  assumptions  are  needed. 

5 . 4 . 1 . 2 Elimination  of  the  Effect  of  Aircra ft  Moti on  from 
the  Doppler  Shifts.  The  first  assumption  deals  with  how  the 
effect  of  aircraft  motion  is  eliminated  from  the  Doppler 
shifts.  Equation  (5. IS)  may  be  expressed  as 

R = (radial  component  of  V)  - VA  cos  (V'A,R) 

(5.20) 


a v v 


-•  v.  v. 
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where  (V^,R)  denotes  the  angle  between  the  aircraft  velocity 
vector  and  the  range  vector  to  the  target.  The  radial  com- 
ponent of  V,  the  target  velocity,  is  what  is  wanted;  the 
other  term,  which  represents  the  effect  of  the  aircraft 
motion,  is  an  error  term  to  be  eliminated  to  the  extent  pos- 
sible. The  extent  to  which  the  elimination  of  the  error 
term  is  possible  is  the  subject  of  the  first  assumption.  The 
way  in  which  the  error- term  reduction  is  performed  may  be 
considered  to  be  by  control  of  the  receiver  local-oscillator 
frequency.  To  eliminate  the  error-term  completely  would  be 
possible  if  only  a single  point  were  illuminated  at  any 
instant,  giving  a unique  value  (over  the  resolution  cell) 
for  the  correction  term  V^cos  (V^,R).  Even  in  that  case, 
the  local-oscillator  frequency  would  have  to  be  modulated 
rapidly.  It  is  assumed  that  this  rapid  modulation  is  not 
feasible.  The  angle  (V^,R)  varies  with  the  range  sweep, 
but  its  variation  with  the  azimuth  beam  angle  is  much 
larger  and  much  slower.  It  is  assumed  that  the  correction 
will  be  synchronized  with  the  scan  of  the  antenna  but  not 
with  the  range  sweeps.  At  each  beam  position,  therefore, 
the  correction  will  be  exact  for  one  target  position,  taken 
to  be  on  the  ground  at  the  middle  of  the  azimuth  beam  angle 
halfway  between  minimum  and  maximum  range.  The  resulting 
error  in  Doppler  shift  is  then  determined  for  any  other  posi- 
tion of  the  echoing  object. 

5.4. 1.3  Nature  of  the  MTI  Filter.  The  second  assumption 
deals  with  the  nature  of  the  MTI  filter.  It  is  assumed  that 
the  filter  is  of  the  form  described  in  Modern  Radar*  as 


* Bcrkowitz , R.  S.  Modern  Radar:  Analys i s , Evaluation,  and 
System  Design.  John  ikilev  f,  Son,  New  York)  (1965)  . 


composed  of  single-delay  filters  in  cascade.  The  filter  out- 
put depends  on  the  filter  transfer  function  as  well  as  on  the 
spectrum  of  the  input.  The  second  assumption  is  dictated 
by  a desire  to  keep  the  program  options  manageable  and  yet  to 
permit  some  choice;  in  this  case,  the  number  of  delay  lines. 
The  transfer  function  is  then  the  function  of  a single 
parameter  N,  the  number  of  delay  lines. 


5.4.2 


Detailed  Computations 


This  paragraph  contains  the  equations  for  computing 
(1)  the  spectrum  of  the  received  echo;  (2)  signal  power  (S) , 
clutter  power  (C) , and  noise  power  (N)  input  to  the  sensor 
and  (3)  S,C,  and  N output  by  the  sensor. 


5. 4. 2.1  Frequency  Spectrum.  The  resolution  cell  extends 
in  the  X- direct ion  over 


RB/2  < XT  < XA  + RB/2 


(5.21) 


where  the  subscripts  T and  A refer  to  the  target  and 
aircraft,  respectively,  R is  the  slant  range  to  the  target, 
and  B is  the  beam  width  angle  in  the  radians. 


Doppler  frequency  runs 


where 


-f  < f < f 

(5.22) 

m - — m 

f = Bf  V.  / c 
m r A 

(5.23) 

is 

aircraft  velocity 

is 

the 

speed  of  light 

is 

the 

radar  tiansmitted  frequency 

W V. 


For  each  target  there  must  be  computed  the  target 
Doppler  frequency  for  entering  the  filter.  The  radial 
velocity  component  is: 


vTR  = (i/R)  [ cxT  - xA)  (vTX  - vA)  + (yt  - ya)  vty 


(5.24) 

(assuming  = 0)  where  - X^  is  the  instantaneous 

distance  component  along  line  of  flight,  - Y^  the  cross- 

flightline distance  component,  R the  slant  range,  V^x 
and  V,py  the  target  velocity  components,  and  VA  is  aircraft 
ground  velocity.  Then  for  the  Doppler  frequency  shift  enter- 
ing the  filter,  associated  with  target  return, 


fT  = -2  (fu  / c)  V„ 


(5.25) 


If  there  is  precipitation,  the  resolution  cell 
volume  will  also  produce  an  echo  from  the  precipitation, 
with  a continuous  power- dens ity  versus  frequency  character- 
istic. Three  frequencies  are  calculated  because  the  total 
clutter  power  is  partitioned  into  three  parts,  each  associated 
with  a vertical  beam  angle  segment  4>j  5 then 


= arcs  in  [ (H.  - HT)  / R ] 


(5.26) 


where  <J>  is  the  depression  angle  for  the  top  of  the  radar 


beam,  and 


<f>.  = <{>  . + (l)  (4>,  - 4>  . ) , i = 0,1,2 

j mm  3 nun'  • 


(5.27) 


For  j = 0,1, 2, 3 set: 


X.  - X. 
J A 


(XT  - XA)  R cos  <}). 


[(XT  - xA)2  ♦ (YT  - ya)2]  3 


(5.2S) 
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and 


Y.  - Y.  = 
J A 


(Yt  - Ya)  R cos  d>3. 


[CXT  - xA);  * (yt  - ya)!] 


2T  -S 


(5.29) 


and 


Z . = R sin  4> . 
J J 


(5.30) 


Let  Yp^,  Vpy,  Vpz  be  the  velocity  components 
of  the  precipitation  VpZ  positive  downward.  Where  Vp^ 
and  Vpy  are  the  wind- velocity  components  and  VpZ  is 
one  input  constant  for  rain  and  another  for  snow.  Let: 


vp.j  ■ &[cw  vpx  * Cy3-YA)  vpy  * Vpz. 


and 


VPE , j = VA  (Xj  - V WR  ' VRC) 


(5.31) 

(5.32) 


where 


1/RC  * 0.5  (1/R  + 1/R  . ) 

v max  rniir 


(5 . 33) 


Finally,  the  three  frequencies  associated  with  the  partitioned 
clutter  power  are: 


fPC,j  ~ r(fr/<0  ('pj  + VPE,j  + VP,j-l  + VPE , j - 1^  ’ ;i=1,2,3 

(5.34) 


5. 4. 2. 2 Power  In.  All  power  return  can  be  characterized  by 
the  general  equation: 

PR  = UV  (5.35) 

where  V is  characteristic  of  the  source  of  the  return  and 
U is  common,  regardless  of  the  source,  and  is  given  by: 


U 


PTpCRn?X?Bc 


4ttB2  (cot  (o^  - cot  uo) 2 (H^-II.j.)  4 


(5.36) 


where 


Pp  = peak  pul  e power 

^CR  = Pu^sC  c<-  pi'ession  ratio 
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n = antenna  efficiency 

A = radar  wave  length 

e = precipitation  attenuation 

B = horizontal  beam  width 

= maximum  depression  angle  = 90° 

w0  = minimum  depression  angle 

= aircraft  height  above  mean  ground  level 

H,p  = target  height  above  mean  ground  level 

3 = the  effective  enhancement  due  to 

sweep-to-sweep  processing.  3 is  found 
by  letting  N be  the  effective  number 
of  sweeps  processed.  (If  uniform 
weighting  is  employed  N is  actual. 

If  shaded  weighting  is  employed  (e.g., 
linear)  N is  about  half  of  actual, 
if  parabolic  weighting,  N is  two- 
thirds  of  actual.  However,  for  input 
we  take  the  effective  N directly.) 

For  incoherent  processing  including 
normal  MTI 

B = /~N 

For  coherent  processing 
3 = N 

For  pseudo- coherent  processing  (i.e., 
video- coherent  or  phase-insensitive 
processing)  the  /— r-?  case  should  be 
used . 


To  compute  V for  each  case,  the  factors  AR  and 
A^,  arc  needed.  AR  is  the  ground  illuminated  area  for  a 
background  cell  and  A,p  corresponds  to  a target  cell. 

AR  = BR  [ (p2  + h (2R+h) ) ^ - p]  (5.37 


where 


i 


l 


B = horizontal  beam  width 
R = slant  range  to  target 
p = ground  range  to  target, 
h = radar  range  resolution 


. 

V 

V 


S 
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At  = aT  / rT 


where 


o^,  = the  target  cross  section 
rT  = target  reflectivity 

The  various  power  contributions  may  now  be  computed.  Power 
corresponding  to  receipt  of  target  return  in  a target  cell 
is : 

PTT  = Ur^  Min  { AR,  A^  } watts  (5.39 

Power  corresponding  to  receipt  of  background  return  in  the 
target  cell  is: 


= UrR  (Ar  - Min  {AR,A.p})  watts 


(5.40 


where  rg  is  the  reflectivity  of  the  background.  The 
background  resolution  cells  will  produce  in  the  area  of  a 
target : 


Pb  - ^Tg  Ar  watts 


(5.41 


5.4. 2. 3 Power  Out.  The  power  coming  out  of  the  filter  for 
each  component  of  power  returned  to  the  resolution  cells  may 
be  computed,  so  that  the  power  spectrum  C(f)  for  the 
clutter  from  the  background  is  as  shown: 


Figure  5-5 

Power  Spectrum  - Clutter  From  Background 
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The  clutter  power  coining  out  of  the  MTI  filter  is 


QB  = | C (f ) G (f ) df  = (PB  / 2fm)  p G(f)df  (5.4 


Where  the  power  transfer  function  G(f)  is  given  by 


G(f)  = (4  sin2  irfx) 


(5.4 


where  N is  the  number  of  delay  lines  and  t is  the  pulse 
repetition  period. 


G (f ) df 


dN  r . 

= 4 J si 


m sin2N  (Tvxf)df  = ^ [ msin2Nf  df 

ITT  j 


-ITT  f 


and,  for  nrf  <<  1 as  is  always  the  case,  sin  f s f and 


fm  4N  2(7iTf  )2N+1 

G(f)df  = ip S 

J 2N  + 1 


(5.4 


then 


Qn  " P„  watts 

2N  + 1 h 


The  filter  output  signal  power  is  for  a resolution  cell 
containing  the  target  at  XT,  YT: 


f 2 it  r f r ’ 

Qt  = PTT  G (fT)  + — — — PT  watts  (5.4 

2N  +1  11 


Signal  power  resulting  from  precipitation  clutter  is: 


Qp  ^PCl^^PCl^  + + (5.4 


PC2  ^ PC2'  PC3  1 PC5J 
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In  equation  (.5.46)  of  the  cited  AIRS  reports,  CG  is  to  be 
identified  with  Pg  above  in  equation  (5.26)  for  present 
purposes.  Then 

PpCj  “ ^2  E(c°t  (2  + esc 

- (cot  4-j)  (2  + CSC2  ^j)]  , j = 1,  2,  3 

(5.48) 


where 


K2  “ PB  aR 


H4  / 3 ar  R' 


(5.49) 


using  the  symbols  of  (3.46).  Th 
clutter  power  into  three  parts 


is  partitions  the  total 


P = P + P + P 
PC  PCI  PC2  PC3 

ich  associated  with  a vertical  beam  angle  segment. 


(5.50) 


m 


To  compute  noise  power  out  (Q^)  , a noise  power 

input  is  required  in  watts  referred  to  the  same  point  as  the 

powers  P of  equations  (5.39),  (5.40),  and  (5.48).  This 

power  is  called  P^.,  in  watts,  spread  evenly  over  the  system 

bandwidth:  „ , . „ ..... 

PN  = (4  E- 15)  (Bl\ ) (xNF)  watts  (5.51) 

where 


Then 


Bh  = bandwidth  in  N!II 

z 

XF  = noise  figure,  absolute  units 


Qn  = PN  (average  value  of  G(f))  = P 


(5.52) 
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5. 4. 2. 4 Ratio  of  Powers.  Combining  from  above  for  signal 
power 

s = QT  + QN  + QP 

for  clutter-plus-noise  power 

(ON)  = QB  + QN  + Qp 

Final ly : 

S/(C  + N)  = (QT+QN+QP)  / (Qb+W  (5-55) 

is  the  argument  of  the  cubic  polynomial  of  (3.103)  of  the 
cited  AIRS  report. 

5 . 5 MOVING  TARGET  DETECTION  BY  FORWARD-LOOKING  RADAR 

5.5.1  Introduction 

The  detection  of  moving  targets  by  forward-looking 
radar  model  (MTIFLR)  has  been  refined  by  the  more  precise 
calculation  of  previously  approximated  equipment  and  environ- 
mental variables.  The  calculations  permit  the  model  to  take 
into  account  the  effects  of  aircraft  speed  and  precipitation 
particle  velocity.  To  achieve  these  refinements  the  user 
must  provide  about  twice  as  many  input  parameters  as  were 
needed  for  the  simpler  model. 

Paragraph  5.5.2  contains  the  equations  for  the 
detailed  computations.  Paragraph  5.5.3  discusses  the  pro- 
gramming changes  to  the  existing  AIRS  SCENARIO  and  EXECUTIVE 
computer  programs.  The  additional  input  parameters  required 
by  the  changes  to  the  model  are  given  in  Paragraph  5.5.4. 
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5. 5. 1.1  Outline  of  Revised  Computations.  The  revised 
general  computations  for  MTIFLR  are  identical  to  the  changes 
for  MTISLR  discussed  in  Paragraph  5.4.1.  Recause  of  dif- 
ferences in  the  geometry  of  MTIFLR  and  MTISLR,  there  arc 
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several  differences  in  the  detailed  computations.  The 
equations  for  the  detailed  computations  for  MTIFLR  arc 
presented  in  Paragraph  5.5.2.  Many  paragraphs  are  the  same 
for  the  two  sensors,  but  the  authors  judge  the  redundancy  to 
be  preferable  to  extensive  cross  references. 


5.5.2 


Detailed  Computations 


This  paragraph  contains  the  equations  for  computing 
the  spectrum  of  the  received  echo  (Paragraph  5. 5. 2.1);  signal 
power  (S) , clutter  (C) , and  noise  power  (N)  input  to  the 
filter  (Paragraph  5. 5. 2. 3). 


5. 5. 2.1  Frequency  Spectrum.  For  each  target  the  target 
Doppler  frequency  for  entering  the  filter  must  be  computed. 
The  radial  velocity*  component  is: 


vTr  C'R-)rCXT*XA-)  ^VTX‘VA^  + <YT-YA>  VTY^ 


(5.56) 


To  V,j.R  must  be  added  the  error  in  correcting  for 

aircraft  motion. 


where 


VTRE  • VW  O/R-l/Kc) 


1/R„  = 1/R  + 1/R  . 

C max  min 


(5.57) 


(5.58) 


Then  for  the  Doppler  frequency  shift  entering  the  filter, 
associated  with  target  return, 


fT  = " 2 (f  / C)  (VTp  + Vrppp) 


(5.59) 


Assuming  V^R  = 0,  Xy-X^  *s  t^e  in?tantancous  distance  com- 
ponent along  line  of  flight,  Y.^-Y^  the  cross  - fl  i ghtline 
distance  component,  R the  slant  range,  and  V^y 

the  target  velocity  components,  and  is  aircraft  ground 

velocity.  R and  R . arc  radar  range  limits. 

' max  nun  ^ 
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where 


R 


the  speed  of  light 

the  radar  transmitted  frequency 


The  frequency  associated  with  background  return  in  the 


resolution  cell  is  fg: 


fB  = -2(fr/c)  V, 


TRE 


(5.60) 


If  there  is  precipitation,  the  resolution  cell 
volume  will  also  produce  an  echo  from  the  precipitation, 
with  a continuous  power-density  versus  frequency  character- 
istic. Three  frequencies  are  calculated  because  the  total 
clutter  power  is  partitioned  into  three  parts,  each  associa- 
ted with  a vertical  beam  angle  segment  <j >y  Then: 


4)^  = arcsin[(H.-HT)  / R] 


(5.61) 


and 


♦j  • Vn  * 'VW-  i ' °>1-2  <s-62’ 


where  is  the  depression  angle  for  the  top  of  the  radar 

beam. 


Set,  for  j = 0,  1,  2,  3 


VXA  " tw  R cos  qftw2  * (5-63) 

Yj-YA  - R cos  V(vv' 4 'wt-  (S>#4) 

Zj  = R sin  4 j (5.65) 


and 


and 
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Let  Vpj,,  Vp^,  Vp ^ be  the  velocity  components  of 
the  precipitation  (Vpz  positive  downward).  Where  Vpx  and 
Vpy  are  the  wind- velocity  components  and  Vpz  is  one  input 
constant  for  rain  and  another  for  snow.  Then,  let: 

vpj  ■ &ftyy  vpx  - tyv  vpy  * ypzi  ^5-66) 


vpej  - vyv  o/r  - i/rc) 


(5.67) 


The  three  frequencies  associated  with  the  partitioned  clutter 
power  are: 

fPC,j  = " (fr/c) (Vp, j+VPE, j+VP,j-l+VPEf j-1^  j =1,2, 3 


(5.68) 


5. 5. 2. 2 Power  In.  All  power  return  to  the  filter  can  be 
characterized  by  the  general  equation: 


v-.l 


PR  = UV 


(5.69) 


wliere  V is  characteristic  of  the  source  of  the  return  and 
U is  common,  regardless  of  the  source,  and  is  given  by: 


where 


PTPCRn2*?ee 


4ttB2  (cot  0)^  - cot  u0) 2 (ha  ' H-p) 


(5.70) 


Pp,  = peak  pulse  power 

PCR=  pulse  compression  ratio 

n = antenna  efficiency 

X = radar  wave  length 

e = precipitation  attenuation 

B = horizontal  beam  width 

Up  = maximum  depression  angle  = 90° 

u0  = minimum  depression  angle 

lip,  = target  height  above  mean  ground  level 
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= aircraft  height  above  mean  ground  level 

B = the  effective  enhancement  due  to 

sweep- to- sweep  processing.  B is  found 
by  letting  N be  the  effective  number 
of  sweeps  processed.  (If  uniform 
weighting  is  employed,  N is  actual. 

If  shaded  weighting  is  employed  (c.g., 
linear),  N is  about  half  of  actual. 

If  parabolic  weighting,  N is  two- 
thirds  of  actual.  However,  for  input 
we  take  the  effective  N directly.) 

For  incoherent  processing, including 
normal  MTI , 

B = /TT 

For  coherent  processing 
B = N 

For  pseudo- coherent  processing  (i.e., 
video-coherent  or  phase  - insens i tive 
processing)  the  /-r.-  case  should  be 
used . 


In  order  to  compute  V for  each  case  the  factors 
AR  and  A,^  are  needed.  A{)  is  the  ground  illuminated  area 
for  a background  cell  and  A^  corresponds  to  a target  cell. 

Ar  = BR  [ (p2  + h (2R+h))!i  - p]  (5.71) 

where 

B = horizontal  beam  width 
R = slant  range  to  target 
p = ground  range  to  target 
h - radar  range  resolution 


where 

Aj  = 

0^  /Tj 

°T 

the  target  cross  section 

rT 

target  reflectivity 

(5.72) 


Finally,  the  various  power  contributions  may  now  be  computed. 
Power  corresponding  to  receipt  of  target  return  in  a target 
cell  is: 


(5.73) 


WJV  J^vr»  VL-.I.  » ff'.'-'-^ 


% L">  k “ k ■ U - ■ 


P tt  = UrT  Min  { AR,  A^.  } watts 


Power  corresponding  to  receipt  of  background  return  in  the 
target  cell  is: 


PTB  = UrB  (Ar  - Min  {AR,AT>)  watts 


(5.74) 


B is  the  reflectivity  of  the  background.  The 


where  r^ 

background  resolution  cells  will  produce  in  the  area  of  a 
target: 


P B = UrRAR  watts 


(5.75) 


5. 5. 2. 3 Power  Out.  The  power  coming  out  of  the  filter  for 
each  component  of  power  returned  to  the  resolution  cells  may 
be  computed.  The  filter  output  signal  power  is  for  a 
resolution  cell  containing  the  target  at 


PTT  G(fT)  + PTU  G ( f R ) 


TB 


B- 


(5.76) 


where  the  power  transfer  function  G(f)  is  given  by 


G (f ) = (4  sin2  Trf t) 


N 


(5.77) 


where  N is  the  number  of  delay  lines  and  t is  the  pulse 
repetition  period. 


The  filter  output  corresponding  to  background 


return  is: 


Qb  = PR  G(f^)  watts 


(5.78) 


Signal  power  resulting  from  precipitation  clutter  is: 


Qn  - Ppn  G (fPn  ) + ^nr?G  (f  ) + Pnp-zG  (fnr,) 


(5.79) 


In  equation  (3.46)  of  the  AIRS  reports,  is  to  be  identi 


fied  with  P above  in  equation  (5.75),  for  present  purposes. 
B 
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Then 


PPC-  = K2  ^cot  *j-i)  + CSc2  *3-1) 

- (cot  ) (2  + esc2  4>^. )]  , 

V 

where 

K2  PB  aR  H / 3 aG  R 

using  the  symbols  of  (3.46).  This  partitions  the 
clutter  power  into  three  parts 

P = P + P + P 
*PC  *PC1  1 PC2  *PC3 


(5.80) 

(5.81) 


total 


(5.82) 


each  associated  with  a vertical  beam  angle  segment. 

To  compute  the  noise  power  out  (Q^) , a noise  power 
input  in  watts  referred  to  the  same  point  as  the  powers  P 
of  (5.73),  (5.74),  and  (5.75)  is  needed.  This  power  is 
called  P^,  in  watts,  spread  evenly  over  the  system  band- 
width 

PN  = (4  E- 15)  (BW)  (NF)  watts  (5.83) 


where 

BW  = bandwidth  in  Mil 

z 

NF  = noise  figure,  absolute  units 


then 


Average  value  of  G(f)) 


(2N)  ! p 

:)2 


(5.84) 


5. 5. 2. 4 Ratio  of  Powers.  Combining  from  above  for  signal 
power 

S ■ Qt  + Qn  + Qp  (5.85) 

for  cluttcr-plus-noise  power 


(C  + N)  - Qn  + Qm  + QP 


(5.86) 


r 


i 

< 
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Finally: 

s/cc  + x)  = (QT+QN+QP)  / (QB+QN+QP)  (5.87) 

is  the  argument  of  the  cubic  polynomial  of  (3.103)  of  the 
cited  AIRS  report. 


5.6  ECM 

5.6.1  Introduction 

The  ECM  model  has  been  extended  to  include  the 
possibility  that  a signal  .cannot  be  detected  due  to  the 
density  of  the  EM  environment. 

Though  not  implemented,  one  alternative  approach 
to  how  this  could  be  modeled  was  presented  in  the  AIRS  Report 
Vol.  I,  Paragraph  3. 2. 5.1).  This  alternative  was  rather 
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5.6.2 


Determination  of  Average -EM  Environment 


Characteristics 
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The  statistics  of  the  average  EM  environment  are 
computed  from  data  on  each  of  the  radars  in  the  scenario 
being  studied.  In  particular,  the  following  statistics 
represent  an  average  EM  environment: 

2 

The  mean  radar  density  per  NM 

G^j  The  mean  radar  transmitter  gain  in  db 
The  mean  PRF  of  the  radar 

H The  ECM  horizon  in  NM 

X The  expected  number  of  interfering 
radars  at  any  time 

a The  expected  value  of  the  fractional 


Where : 


i is  an  index  over  all  radar  types 

n.  is  the  number  of  radars  of  type  i in 
1 the  scenario 

w.  is  the  percent  of  time  a radar  of 
type  i emits 

g.  is  the  transmitter  gain  for  radar 
type  i 

r^  is  the  PRF  of  a radar  of  type  i 

in  order  to  find  the  mean  radar  density,  we  must 
first  find  the  area,  denoted  by  A,  in  NM  covered  by 
the  relevant  scenario.  It  will  be  recalled  (see  page  2-30 
of  AIRS,  Vol.  I)  that  the  user  provides  the  simulation 
with  the  maximum  and  minimum  longitude  and  latitude; 

let : 

L°MAX  = longitude  in  degrees 

L°min  = Minimum  longitude  in  degrees 

^aMAX  = Maxi-mum  latitude  in  degrees 

L^MIN  = Minimum  latitude  in  degrees 

The  area  determined  by  these  points  is  given  by  the  following 
two  equations: 


LoMAX  " Lo'!IN 


2irr?  cos  0 d0 


(5.89) 
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(Jj0MAX  ' LoMIN)27rr‘ 


^LaMAX^  ’ Sin  ^LaMIN' 


(5.90) 


and  r is  taken  as  [21600/2tt]  nautical  miles 


Thus,  the  density  of  interferring  radars  per  NM  is  just: 


(5.91) 


where : 


n.  is  the  number  of  radars  of  type  i 
in  the  scenario 

i is  the  index  over  all  radar  types 


The  ECM  horizon  is  defined  as  that  distance  at 
which  the  probability  that  a radar  (at  elevation  +2o 
above  mean  ground  level  --  a typical  height  for  radars) 
is  obscured  by  terrain  is  just  0.5.  Paragraph  2.4  of 
AIRS,  Vol.  I,  describes  the  terrain  algorithm  built  into 
the  SCENARIO  model.  As  discussed  in  the  referenced  para- 
graph, the  algorithm  was  run  off  line  and  the  resulting 
curves  incorporated  directly  into  SCENARIO.  The  algorithm 
was  fully  normalized  so  that  the  resulting  probability  of 
non-blocking  (as  a function  of  depression  angle  from  the  tar- 
get) was  not  a function  of  the  aircraft  height,  wavelength  of 
terrain,  or  standard  deviation  of  the  height  of  the  terrain 
about  its  mean.  Using  the  model,  the  following  graph  was 

( 1 n\rn  1 nnrwl  • 


,'vvvy. 


Thus,  for  a target  at  height  +2a  above  mean  ground  level 
(AMGL) , the  probability  that  the  target  is  blocked  is  0.5 
at  a depression  angle  of  0.13  radians.  Therefore,  II  is 
found  by : 


(ha  - 2o) 


ii-  a 7.69 

H - (h  - 2a) 

6080  tan ( . 13)  6080.  a 


(5.92) 


where 


h = height  of  aircraft  above  mean  ground 
a level  in  feet 

a = standard  deviation  of  terrain  height 
about  its  mean 

and  H has  the  units  of  NM. 

The  expected  number  of  interfering  radars  in  view  at 
at  any  one  time  is  given  by: 

X = ttH2Dm 

For  the  average  radar,  the  expected  value  of  fractional 
main  beam  time  is: 


it  10 


'cyioy 


(5.93) 


All  these  parameters  of  the  average  EM  environment  are  com- 
puted in  the  SCENARIO  program.  Of  them,  a,  R^,  D^,  and  X 
are  transferred  to  the  EXECUTIVE  program  via  SCENARIO  output 
tape . 
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5.6.3 


Overload  Analysis 


5.6. 3.1  Assumptions  and  Definitions  of  Terms.  The  follow- 
ing terms  will  be  used  in  the  overload  analysis: 

Sweep  --  the  interpulse  period  of  a radar 

Radar  Scan  --  inter  main  lobe  period  of  a radar 

ECM  Scan  --  one  cycle  of  the  ECM  sweep  receiver 

ECM  Look  --  the  time  period  during  which  a 

radar  is  within  a given  ECM  beam. 

ECM  Dwell  Time  --  the  time  during  which  an  ECM 
receiver  dwells  on  a frequency. 

The  following  assumptions  are  made: 

(1)  The  radar  scan  period  is  much  longer  than 
the  ECM  scan  period. 

(2)  The  ECM  look  period  is  much  longer  than 
the  ECM  scan  period. 

(3)  The  primary  interference  effect  is  due  to 
overloading  of  the  ECM  processor,  which,  upon 
overload  takes  the  first  pul ses presented  to 
it  per  unit  time  provided  they  are  not  so 
different  in  amplitude  as  given  by  the  front 
to  back  ratio  of  each. 

(4)  Mon-marginal  action  - at  a given  instant  the 
radar  signal  is  either  above  or  below  the 
interference  level;  probabilistic  sweep 
chopouts  may  be  ignored. 

Let  the  following  parameters  be  given  for  the  ECM  (via  user 
input  to  EXECUTIVE): 

R,  = maximum  ECM  signal  processing  rate  in 
^ pulses  per  second 

= the  number  of  ECM  receive  channels 

B.  = the  ECM  instantaneous  receive  bandwidth, 

1 in  MU 

z 

B0  = the  ECM  total  system  bandwidth,  in  MH 

L Z 

Finally,  a given  radar  (radar  j)  lias  the  following  parameters: 

T.  = transmitter  gain  in  db 

F.  = ERF 
J 

P„  . = fraction  of  time  emitting  * 

RAD , j h A 
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5. 6. 3. 2 Computation  of  Overload  Probabilities.  R is  the 

p 

maximum  ECM  signal  processing  rate.  Then,  if 

(V  « < RP 

no  interference  will  exist  in  either  the  main  or  side  lobes. 
If  this  condition  does  not  exist,  interference  may  occur. 

Let  R^  be  the  estimated  fractional  dwell  time  of  the  ECM; 
is  given  by: 


R.  = (B.)  (Nd)  / B, 


(5.94) 


Then  Rq  ^ , the  number  of  competing  pulses  per  second  in 
the  main  lobe  for  target  j is: 

V;  ■ Cfj  + © tRM)m  ] Rd  (5.93) 

and  the  non-interference  probability  for  the  main  lobe, 

P is* 

ninta, j ' * 


P . . 
nmta,  j 


if  Rx,X  < R 
M p 


MIN  1 if  RMX  >•  R 

R I>  — p 

a v 


(5.96) 


in  the  side  lobe  for  target  j is  given  by: 

RM  ■ * d-^RyX]  Rd 

and  the  side  lobe  non-interference  probability,  P . 

r ; ’ ni 

/ 1 . 0 if  R..X  < R 


P . _ . 
n i n t B , j 


if  R..X  < R 
M p 


MIN  -rt1-,  1 if  R..X  > R 
]V  M ~ P 


(5.9$) 


second 

yl*-.1 
*.  % 1 
aA 

(5.97) 

v/v 

. bv : 

» J 

'Sv« 

>v 

/ vli 

5-39 


-J  •J.vvv.w.v.v.-.  .V  .v  J 


A 

LX  [X 


v*: 

• 


>> 


rr T7  *■'  •'  v .■  . .•  v 


5.6.4  Computation  of  Target  Detectability 

Additional  information  must  be  computed  at  this 
point,  namely,  the  fractional  main  beam  time  for  target  j. 
Let  this  be  denoted  as  oc^  ; then 


‘j  " 7 'Ij/lOy 
TT  10 


(5.99) 


Nov,  for  each  look  of  the  ECM,  the  detectability  of  target  j, 
denoted  by  PIE  , is  found  by: 

PD.  = PMC  . PD.n  . fa-PDC  .P  . . - + (1 - a • ) PDC„  .P  . .a  .] 

3 NS, 3 PAD, 3 l 3 a, 3 ninta.3  3 3 , J mnt6,3j 


where 


RAD,  3 


NS, 3 


(5.100) 

= Pr  (target  3*  is  emitting),  i.e.,  percent  of  time 
target  3 is  emitting. 

= the  "raw"  detectability  [i.e.,  assuming  (1)  the 
target  was  emitting,  (2)  the  target  was  not  ob- 
scured by  terrain,  and  (3)  the  target  was  not 
lost  due  to  overloading]  in  the  main  lobe  for 
target  3 . 

= the  "raw"  detectability  of  the  side  lobe  for 
target  3 . 

= the  probability  that  target  3 is  not  obscured 
on  the  present  look  by  terrain. 


PDC  . and  PDC„  . are  computed  as  shown  in  AIRS,  Vol.  I. 

a,  3 p > J rp-j 

Paragraph  3. 2. 6.1  as  P^  and  respectively.  Let  PD. 

(k  = 1,  2,  ...)  be  the  detectability  of  target  3 on  look  k. 
This  set  is  found  via  equation  (5.100)  computed  once  for  each 
look.  Then  the  overall  target  detectability  for  target  3, 
denoted  by  PD j , is: 


PD.  = fl-n  (1-PD.  (-KJ)]  P 

J J J up 


(5.101) 


where  P is  the  probability  that  the  IiCM  equipment  is 
operational . 
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5.6.5  Identi fiability  and  Loca 1 i zab i 1 i ty 

The  methods  discussed  in  AIRS,  Vol.  I for  computing 
the  identif iabil ity  and  localizability  of  a given  target  are 
noj:  changed. 

5.7  PHOTO 

5.7.1  Introduction 

A number  of  improvements  were  made  in  the  PHOTO 
model  as  described  in  Paragraph  3.2.7  of  AIRS,  Vol.  1. 
Specifically,  the  methods  for  computing  contrast,  for  both 
day  and  night  photography,  have  been  improved.  Paragraph 

5.7.2  of  this  report  gives  the  new  method  for  computing  con- 
trast for  daylight  photography , while  Paragraph  5.7.3  pre- 
sents contrast  computation  for  night  photography. 

5.7.2  Contrast. -Day  Photography 

The  method  given  in  Paragraph  3. 2. 7. 1.1  of  AIRS, 
Vol.  I,  for  computing  the  logarithmic  contrast  has  been 
modified.  Let: 

r = target  reflectivity 

rK  = background  reflectivity 

u - 1 
X^  = scattering  coefficient  for  haze  --  NM 

X = scattering  coefficient  for  rain  --  NM  1 

r - 1 
Xa  = scattering  coefficient  for  clear  air  --  NM 

R^  = slant  range  to  the  target  through  haze  --  NM 

Rr  = slant  range  to  the  target  through  rain  --  NM 

R = slant  range  to  the  target  through  clear 

air  --  NM 

Bo  = brightness  of  the  sky 
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Consider  a target  of  reflectivity  r^ , illuminated  by  a sky 
of  uniform  brightness  B . The  target  is  considered  small 
(that  is,  it  is  assumed  to  subtend  less  than  5°  of  arc)  so 
that  second  order  scattering  may  be  ignored.  Assuming  opera- 
tion on  the  linear  portion  of  the  film  gamma  scale,  the 
perceived  brightness  of  the  target,  denoted  by  B^,  is  given 
by: 


B . = B 
t o 


t -(WWW 


(5.102) 

In  this  equation,  the  first  term  on  the  right  of  the  equality 
represents  scattering;  the  second  term  represents  attenuation. 
Absorption  is  ignored.  Similarly,  the  preceived  brightness 
of  the  background,  denoted  by  B^,  is: 


B,  = B 
b o 


n -(WWW,  t -(WWW 

(1-e  ) + rbe 


(5.103) 

In  the  above  equations  for  the  brightness  of  the  target  and 
background,  A'  can  be  taken  as  zero  since  optical  scattering 

3. 

through  clear  air  is  negligible. 

Therefore,  the  logarithmic  contrast,  C is  given  by: 
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C = Y log 


= Y 


10  B, 


1 


log. 


(l-rt), 
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where  y is  the  film  gamma. 


(5.10-1) 


%'V 


• . 
tv*. 


*►  .> 


One  furtlicr  detail  remains;  A.  and  A^  must  be 
developed.  Both  arc  functions  of  the  prevailing  visibility 
(haze  or  rain  as  appropriate). 
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The  following  relationship  is  used*: 

{h  for  haze 
r for  ram 

is  the  scattering  coefficient  and 


(5.105) 


where 


A. 

1 


V. 

1 


is  the  pre- 
vailing visibility  in  statute  miles  (A^  has  the  units  of  NM  1 ) 


Note  that,  in  the  above  development,  shadow  contrast 
has  been  ignored  and  that  scattering  is  assumed  isotropic. 

The  former  assumption  is  good  on  a hazy  day  or  under  clouds; 
it  may  be  taken  to  represent  a worst  case.  Optical  haze 
scattering  is  actually  anisotropic,  with  the  angular  dependence 
highly  sensitive  to  the  distribution  in  size  of  haze  particles. 
Theoretically,  the  appropriate  model  is  that  of  Mie  scattering; 
however,  since  the  particle  size  distribution  is  not  known, 
and  since  experimental  data  yield  variations  of  not  more  than 
0.2  C-units,  these  can  be  ignored. 


5.7.3 


Contrast  - Local  1 1 luininat  ion** 


Consider  the  case  of  loca]  illumination  only  where 
the  ambient  plays  no  role.  The  chief  culprit  of  contrast 
reduction  is  back-scattered  light. 

Terms  to  be  used  in  this  discussion  are  shown  in 
Figure  5-7;  the  light  source  is  assumed  offset  some  distance 

d from  the  line  of  sight  at  an  angle  G from  that  line  of 

sight.  The  angle  0 can  be  taken  as  constant  over  the  field 

of  view  of  the  camera***.  For  simplicity,  d is  assumed 


* Wolfe,  William  L.  (ed.)  Handbook  of  Mi  1 i tary  In  fra red 
Technology , Office  of  Naval  Research,  Dept,  of  theNavy, 


Washington,  D.C.,  1965,  p.  205. 

* * This  model  is  appropriate  for  either  (1)  night  photography, 
or  (2)  monochromatic  photography  with  a laser  source. 

***In  the  implementation,  0 is  taken  to  be  90°. 
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constant  with  time;  such  would  not  be  the  case  if  the  light 
source  were  a flare.  This  assumption  is  tantamount  to  assuming 
the  strobe  lighting  is  used  for  night  photography;  this 
assumption  is  most  likely  correct  with  modern  reconnaissance 
aircraft  where  room  (e.g.,  to  carry  flares)  is  at  a premium. 

As  illustrated  in  Figure  5-5,  the  "field-of-view" 
of  one  resolution  element  is  ft  steradians  (a  typical  value 
of  ft  is  10  ^)  . It  is  further  assumed  that  there  is  a 
minimum  backscattering  range  of  R ; this  minimum  range  would 
be  provided  by  either:  (1)  baffling,  or  (2)  an  electronic 

shuttering  system. 

Let  the  light  source  intensity  in  the  target  direc- 
tion be  L lumens/ster adian  Then  the  illumination  at  the 
target  element  within  ft  (letting  d <<  R) , denoted  by  1^, 
is  given  by: 

I„  = P.R2e"XR  = Lfte'AR  (5.106) 

i 

vheie  ^ slant  range  to  the  target  = R +R^+R 

A = scattering  coefficient 

Ais  taken  as  the  weighted  average  of  A , A^ , and  Af: 


-V 

£ 


A 


RA  + R>  A,  + R A 
a a h h r r 


R 


(5.107) 


The  illumination  received  at  the  camera  from  the  target, 
denoted  by  J.j. , is: 


T t *AR,1  v . ,1s 

J.p  I.pC  ( o)lp0i(  ) 

2 it 


(5.108) 


where  a is  the  camera  aperture  area.  Rewritting: 


T Lft  - 2 AR  ] . 

JT  - — c r a( — ) 

J R“  r 2n 


(5.109) 


5 - -I  5 


vV.-.tV 


r.  v, 


Similarly,  the  illumination  received  at  the  camera  from  the 
background,  denoted  by  J is: 

T _ Ln  - 2 A R , 1 . 

JB  " 72  e rba(~ 3 (5.110; 

K 2 7f 


Backscattering  must  now  be  introduced.  Take  an 

elemental  unit  of  path  (a  thin  cylinder)  at  an  arbitrary 

distance  x and  with  a cross  section  fix2  and  length  dx. 

Its  total  illumination,  I . is: 

P’ 


Ip  = Lfte  cos  $ — 

where 


^(x  - d cos  6) 2 + (d  sin  G)2 
= the  source  - element  - camera  angle  and  is 


arctan 


sin  9 
- cos  0 


This  path  length  scatters  a fraction  Adx  into  a sphere; 
hence,  the  backscattcr  received  by  the  camera,  J.,  is  found 
from:  S 

j t __Lj7  -AfZ+x).  cx 

dJ s - — e **  cos  $ — dx 

2 4 ir 


7 _ LftqA  e ^ 2 x 3 cos  4> 

s , 2 

4 it  I Z 


(5.112; 


This  integral  cannot  be  solved  explicitly,  so  a Lagucrre 
integration  technique  is  used. 


».■» :w  v.VlV.va^v.v  v.w.w  .vjn 


First,  define  a function  G(X)  as: 


Then 


1 _ Lfia  A 

Js  - 
4 7T 


G (X)  = e" (z*x) * cos  q/z2 


G(x)e*2Xx  dx 


G (x) e" 2 Xx  dx 


(5.115) 


Let  y be  defined  as  the  term  in  square  brackets  above; 
further,  let: 


F(k)  = G(x)e'2Xx  dx 


(5.114) 


This  function,  F(k),  is  evaluated  by  the  method  described 
on  page  923  of  the  NBS  Handbook  of  Mathematical  Functions”; 
namely,  the  Laguerre  numerical  integration: 


,°o  n 

e"xf (x)  dx  * J wA  f(x.) 

* 


(5.115) 


where  (w^,  x^)  are  tabulated  as  a function  of  n. 

Using  this  approximation,  and  transforming  the 


variables : 


F(k)  = e"x2X  G (x)  dx  = e'2Xk  \ G + , 

i = l I 2 A I 2A 


In  the  implementation,  n is  chosen  as  5. 

Table  5-1  gives  (wi,  for  i = l,2,...5. 


(5.116) 


Abrainowitz  and  Stegun  (ed.),  U.  S.  Department  of  Commerce, 
lyindbook  oT  Mathematical  Functions,  l).  S.  Government  Printing 
'-'Tficc,  Washington,  D.C.,  June  1964  , page  923. 
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SECTION  VI 


OUTPUT  OF  FACTORS 
AFFECTING  SENSOR  PERFORMANCE 


6.1  INTRODUCTION 

This  section  discusses  a program,  designated  as 
INTER,  which  has  been  developed  to  enable  the  AIRS  user  to 
study  intermediate  sensor  performance  factors  in  an  organ- 
ized fashion. 

6 • 2 THE  INTER  PROGRAM 

Until  this  program  was  developed,  the  user  could 
only  analyze  sensors  by  their  macro -performance ; that  is, 
he  could  examine  the  detectability,  identifiability  and 
localizability  of  each  target -sensor -pass  combination  (by 


6-1 


, -TV  r ■ *'■  * “ 


use  of  the  EXECUTIVE  hard  copy  output --see  paragraph  3.6.3 
of  AIRS,  Volume  1)  or  he  could  examine  how  a sensor  per- 
formed vis-a-vis  groups  of  targets  (by  using  the  EVAL  out - 
puts--see  paragraph  4.6  of  this  report).  However,  the  user 
was  unable  to  determine  why  the  detectabilities,  identifi- 
abilities,  and  locali tab ilities  were  as  they  were.  Thus, 
for  example,  if  the  IR  sensor  continually  gave  low  detect- 
abilities, the  user  would  be  unable  to  determine  if  it  were 
due  to  poor  target  contrast,  bad  weather,  etc. 

For  each  sensor,  a list  of  intermediate  factors 
which  would  be  most  useful  to  the  user  in  answering  ques- 
tions of  why  sensors  performed  as  they  did  was  prepared- - 
these  lists  are  given  in  Table  6.1.  It  would  have  been 
possible  to  output  these  intermediate  factors  during  the 
running  of  the  EXECUTIVE  program  each  time  a particular 
target-sensor-pass  combination  was  considered.  This  method 
of  providing  the  data  to  the  user  would  have  been  both 
clumsy  and  overwhelming;  rather,  it  was  decided  that  all 
intermediate  factors  would  be  stored  on  the  EXECUTIVE  Target/ 
Sensor  Output  Tape  (see  paragraph  7. 3. 2.1  of  this  report). 
This  tape  would  then  serve  as  input  to  the  INTER  program 
which  would  organize  the  data  and  output  it  in  a succinct 
manner . 

The  next  question  was  how  to  portray  the  data  in  a 
manner  that  would  be  most  useful  to  the  user.  After  consid- 
eration of  a number  of  alternatives,  it  was  decided  that  for 
each  sensor  the  factors  of  interest  would  be  averaged  over 
all  valid*  target-sensor-pass  combinations  and  output  as  a 


* Thus , if  on  a particular  pass  a target  never  foil  within  the 
field  of  view  of  a particular  sensor,  this  would  not  degrade 
the  average  performance  factors. 


ange  interval  in  caouiar  lortu.  me  range  inter 
sensor  are  given  in  Table  6.1.  This  method  of 
reduced  the  data  to  a manageable  amount  while  s 
user  with  sufficient  information  to  enable  him 
sensor  performed  as  it  did. 

ce,  for  each  of  the  14  sensors  (the  dummy  senso 

d) > a table  is  printed  out  by  the  INTER  program 

* 

a sample  output  for  the  FLIR  sensor;  all  othe 
nalogous.  Across  the  horizontal  axis  are  the 
c intervals  for  the  sensor  to  which  the  table 
the  vertical  axis  is  a mnemonic  list  of  the  fac 
Table  6-1  for  that  sensor.  At  the  top  of  the 
mmary  of  the  number  of  target-pass  sightings 
ange  interval  by  target  type.  (The  total  numbe 
s siphtincs  within  a piven  range  interval  is  the 

for  that  range  interval 

for  each  sens 
used  for 


NTER  program  is  th 
XECUTIVE  program, 
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TABLE  6-2 
MNEMONIC  LIST 


SENSOR  MNEMONIC  DEFINITION 

IR  CB1CIR  Contrast  in  Band  1 when  clouds  are  present 

CB1NIR  Contrast  in  Band  2 when  clouds  are  absent 

CB2CIR  Contrast  in  Band  2 when  clouds  are  present 

CB2NIR  Contrast  in  Band  2 when  clouds  are  absent 

CDET1  Conditional  Detectabi 1 ity- - Band  1 

CDET2  Conditional  Detectabi 1 i ty- - Band  2 

CDIR-1  Conditional  Detectability- - level  1 

CDIR-2  Conditional  Detectability- - level  2 

CDIR-3  Conditional  Detectabi I ity- - level  3 

CDIR-4  Conditional  Detectability- - level  4 

TDIR-1  Total  Detectability- -level  1 

TDIR-2  Total  Detectability- - level  2 

TDIR-3  Total  Detectability- - level  3 

TDIR-4  Total  Detectability- - level  4 

TIIR-1  Total  Identi f iabi 1 i ty- - level  1 
TIIR-2  Total  Identif iability- - level  2 

TIIR-3  Total  I dent i f i ability - - level  3 

TIIR-4  Total  Identif iability- -level  4 


SLR  SLRC0N1  Contrast- - level  1 

SLRCON24  Contrast -- levels  2,  3,  and  4 

E0RLEV1  Effective  Output  Resolut ion- - level  1 

EORLEV2-4  Effective  Output  Resolution- -levels  2,  3,  and  4 

ESRLEV1  Effective  System  Resolution- - level  1 

ESRLEV2-4  Effective  System  Rcsol ution- - levels  2,  5,  and  4 

CDSLR-1  Conditional  Detectability- - level  1 

CDSLR-2  Conditional  Detectability- - level  2 

CDSLR-3  Conditional  Detectability- - level  3 

CDSLR-4  Conditional  Detectab i lity- - level  4 

TDSLR-1  Total  Detectability- - level  1 

TDSLR- 2 Total  Detectability- - level  2 

TDSLR-3  Total  Detectabili ty- - level  3 

TDSLR-4  Total  Detec t abi li ty- - level  4 

TISLR-1  Total  Identif iability — level  1 

T1SLR-2  Total  Ident i f iabi 1 i ty- - level  2 

TISLR-3  Total  Ident i fiabili ty — level  3 

TISLR-4  Total  Ident i fiabil ityr- - level  4 
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SENSOR 

MTISLR 


P 


MTIFLR 


MNEMONIC 


PDSPSL 

POSFTSL 


POSFBSL 


PCPWSL 

NPFDLSL 


PRASL 


TNPRSL 

CDMRSL-1 

CDMRSL-2 

CDMRSL-3 

CDMRSL-4 

TDMRSL-1 

TDMRSL-2 

TDMRSL-3 

TDMRSL-4 

TJMRSL-1 

TIMRSL-2 

TIMRSL-3 

TIMRSL-4 


RDSPFL 

POSFTFL 


POSFBFL 


PC  PUT  L 
NPFDLFL 


PRAFL 


F.XNONL 


CDMRFL-I 
CUMRFL-2 
CDMRFL- 3 
CDMRFL-4 


TABLE  6-2 


MNEMONIC  LIST  (CONT.) 


DIFINITION 


Radial  Speed  of  Target  in  knots 

Power  output  (in  watts)  at  the  sensor  due 

to  signal  return  from  the  target 

Power  output  (in  watts)  at  the  sensor  due 

to  signal  return  from  background 

Precipitation  clutter  power  (in  watts) 

Noise  power  (in  watts)  introduced  by  the 
sensor  delay  lines 

Ratio  of  target  signal  plus  noise  plus  precip- 
itation clutter  power  to  background  signal 
plus  noise  plus  precipitation  clutter  power 
Terrain  nonshadowing  probability 
Conditional  Detectability-- level 
Conditional  Detect ability --level 
Conditional  Detectability- -level 
Conditional  Detectability- -level 
Total  Detectability- -level  1 
Total  Detectability- - level  2 
Total  Detectability- -level  3 
Total  Detectabi 1 ity- - level  4 
Total  Identi fiabi lity- - level 
Total  Identi f iab i 1 i ty - - level 
Total  Identifi ability --level 
Total  Identif iability- - level 


1 

2 

3 

4 


j. 

2 

3 

4 


Radial  Speed  of 
Power  output  (in 
to  signal  return 
Power  output  (in 
to  signal  return 


target 


in  knots 
wratts)  at  the  sensor 
from  target 
watts)  at  the  sensor 
from  background 


due 


due 


Precipitation  clutter  power  (in  watts) 
Noise  power  (in  watts)  introduced  by  the 
sensor  delav  lines 


Ratio  of  target 
itation  clutter 
plus  noise  plus 
Expected  number 
p a s s 

Condi tional 
Condi t iona 1 
Condi t i onn 1 
Coiul  i t ional 


signal  plus  noise  plus  precip 
power  to  background  signal 
precipitation  clutter  power 
of  unobstructed  looks  on  this 


Detectabi lity- - level 
Detectabi 1 ity- - level 
De  t c c t ab  i 1 i t y - - 1 e ve  1 
Detect  ability --lev el 
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TABU:  6-2 


MNEMONIC  LIST  (CONT.) 


SENSOR  MNEMONIC  DEFINITION 

y: 

TDMRFL-l  Total  Dot ectab i 1 ity- - level  1 >' 

TDMRFL-2  Total  De tectabi 1 i ty - - 1 evel  2 

TDMRFL-3'  Total  Detectab i 1 ity- - level  3 

TDMRFL-4  Total  Detectability- -level  4 

TIMRFL-1  Total  Identifiabi lity- - level  1 

TIMRFL-2  Total  Identif iability- - level  2 

TIMRFL-3  Total  Identifiability- -level  3 

TIMRFL-4  Total  Identif  iability- - level  4 

\ 

PHOTO  CFOTO  Contrast 

RFOTO-12  Resolution- -levels  1 and  2 

RFOTO-34  Resolution- -levels  3 and  4 

FFOTO  Number  of  frames  on  which  target  appears 

ACDET-12  Conditional  Detectability,  averaged  over  all 
looks- - levels  1 and  2 

ACDET-34  Conditional  Detectability,  averaged  over  all 
looks -- levels  3 and  4 

ATDET-12  Total  Detectability,  averaged  over  all  looks-- 
1 eve Is  1 and  2 

ATDET-34  Total  Detectability,  averaged  over  all  looks- - 
levels  3 and  4 

ATIDEN-12  Total  Identifiability,  averaged  over  all 
looks -- levels  1 and  2 

ATIDEN-34  Total  Identifiability,  averaged  over  all 
looks -- levels  3 and  4 

CDF0T0-1  Overall  Conditional  Detectability- -level  1 

CDF0T0-2  Overall  Conditional  Detectability- -level  2 

CDF0T0-3  Overall  Conditional  Detectability- -level  3 

CDF0T0-4  Overall  Conditional  Det ectab i 1 i ty- - 1 evel  4 

TDF0T0-1  Overall  Total  Detectabili  ty- - level  1 
TDF0T0-2  Overall  Total  Detectabil ity- - level  2 
TDF0T0-3  Overall  Total  Detectabi lity- - level  3 
TDF0T0-4  Overall  Total  Detect abi lity- - level  4 
TIF0T0-1  Overall  Total  Identif iabi lity- - level  1 
T1F0T0-2  Overall  Total  Identifiability- -level  2 
TIF0T0-3  Overall  Total  ldcnt i f iabi 1 i ty- - level  5 
TIF0T0-4  Overall  Total  Identif iability- - level  4 

FLIR  (Identical  with  IR  with  one  addition.) 

SLANT  RAN  Slant  range  (in  NM)  to  target  at  last  valid 
look 
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TABLE  6-2 

MNEMONIC  LIST  (CONT.) 


SENSOR  MNEMONIC 


DEFINITION 


ECM 


SNML 

SNSL 

CDML 

CDSL 

NIPRML 

NIPRSL 

CDECM-I 
CDECM-2 
CDECM-3 
CDECM-4 
TDECM-1 
TDECM-2 
TDECM- 3 
TDECM-4 
TIECM-1 
T I ECM- -2 
TIECM-3 
T I ECM -4 


Signal  to  noise  ratio,  main  lobe  radiation 
Signal  to  noise  ratio,  side  lobe  radiation 
Conditional  Detectability  in  main  lobe 
averaged  over  all  looks* 

Conditional  Detectability  in  side  lobe 
averaged  over  all  looks* 

Probability  of  no  interference  in  main  lobe 
averaged  over  all  looks 

Probability  of  no  interference  in  side  lobe 

averaged  over  all  looks 

Conditional  Detectab i li ty- - level  1 

Conditional  Detectabi 1 ity- - level  2 

Conditional  Detectabi lity- - level  3 

Conditional  Detectabili ty- - level  4 

Total  Detectability- - level  1 

Total  Detectability- -level  2 

Total  Detectability --level  3 

Total  Detectability- -level  4 

Total  Identifiability- - level  1 

Total  Ident i fiabi 1 i ty- - level  2 

Total  Identifiability — level  3 

Total  Identifiability- - level  4 


"£ 

Conditional  here  means  assuming  no  terrain  blocking,  equipr.cn 
failure,  or  interference  from  other  radars. 
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INPUTS  AND  OUTPUTS  OF  THE  SCENARIO  AND 
EXECUTIVE  PROGRAMS 


7.1  INTRODUCTION 

The  purpose  of  this  section  is  to  provide  the  reader 
with  the  complete  input  requirements  of  the  SCENARIO  and 
EXECUTIVE  programs  of  the  AIRS  Model.  In  addition,  this  section 
will  provide  descriptions  of  the  outputs  of  these  programs. 
Paragraph  7.2  will  discuss  the  input  and  output  of  the  SCENARIO 
program;  specifically,  the  input  will  be  discussed  in  paragraph 

7.2.1  and  the  output  in  paragraph  7.2.2.  Similarly,  within 
paragraph  7.3  the  inputs  and  outputs  of  the  EXECUTIVE  program 
will  be  discussed  in  7.3.1  and  7.3.2  respectively. 

7.2  SCENARIO  INPUTS  AND  OUTPUTS 


The  only  required  input  to  the  SCENARIO  program  is 
a user-developed  and  supplied  card  input  deck.  The  output 
consists  of  two  tapes.  These  will  be  discussed  in  the  following 
paragraphs . 
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7.2.1 


Inputs 

The  best  way  to  illustrate  the  form  of  the  new 
input  deck  required  by  SCENARIO  is  by  Table  7.1.  This  table 
serves  a dual  purpose:  1)  it  defines  all  necessary  input 
parameters,  and  2)  it  provides  instruction  for  the  preparation 
of  the  input  deck  itself.  This  table  supercedes  and  replaces 
pages  2-30  through  2-36  inclusive  of  AIRS,  Volume  I. 

7.2.2  Outputs 

SCENARIO  outputs  are  on  two,  or  optionally  three 
magnetic  tapes. 

Output  Tape  1 contains  a descriptive  record  of  target 
sightings  by  each  of  the  15  sensors  that  are  switched  on  for  a 
given  run,  and  serves  as  input  to  the  EXECUTIVE  Model. 

The  second  tape,  Output  Tape  2,  contains  a frame-by- 
frame  descriptive  record  of  what  each  of  the  (up  to  8)  cameras 
on  board  the  aircraft  surveyed.  After  sorting,  this  tape 
serves  as  input  to  the  PI-Q  Model  (See  AIRS  Vol.  I,  Section  V; 

* 

The  optional  tape  , Output  Tape  3,  lists  target  number, 
target  type,  and  the  time  at  which  the  aircraft  flew  past  the 
target  for  each  target  sighting.  This  tape  is  sorted  by  time  and 
then  listed  to  provide  a time-line  record  of  the  mission. 

7.2.2. 1 Output  Tape  1.  The  form  of  Output  Tape  1 has  not  been 
changed  from  that  reported  in  paragraph  2.7.2  of  AIRS,  Volume  I. 
Table  2-1  and  Table  2-2  of  AIRS,  Volume  I are  valid;  however, 
Table  7-2  of  this  report  replaces  Tables  2-5a  and  2-3b  of  the 
cited  AIRS  report. 


vv 


AS 


.v> 

1 


% V 


\ % . 


. - . V 

vVa 


* 

This  tape  is  procuded  if  Sense  Switch  2 on  the  Cl'C  Computer 
Console  is  turned  on. 
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7*2,2'2  Output  Tape  2.  The  second  output  tape  of  the  SCENARIO 
program  has  not  been  modified  in  any  manner.  Hence  paragraph 

2.7.3  of  AIRS  Volume  I remains  completely  valid  and  will  not 
be  restated  here. 

7.2.2.3  Output  lape  5.  As  in  the  case  of  Output  Tape  2, 
this  tape  remains  unchanged.  The  reader  is  referred  to 
paragraph  2.7.4  of  AIRS,  Volume  I,  for  a complete  discussion 
of  Output  Tape  3. 


7.3 

7.3.1 

The  inputs  to  the  EXECUTIVE  program  are  of  two  types: 
1)  tape  and  2)  card.  The  tape  input  to  EXECUTIVE  is  Output 
Tape  #1  of  the  SCENARIO  program,  as  discussed  in  paragraph 

7 . 2 . 2 . 1 above . 

Table  7.3  gives  the  form  of  the  EXECUTIVE  input  deck 
by  card  number,  variable  name,  variable  definition  and  format. 

It  is  sufficiently  complete  to  serve  as  a guide  to  the  AIRS 
user  in  developing  an  EXECUTIVE  input  deck.  There  arc  55  cards 
in  the  input  deck,  each  card(s)  pertains  to  a specific  sensor 
as  follows: 


EXECUTIVE  INPUTS  AND  OUTPUTS 
Inputs 


znx 


AV. : -N  :• 


CARDS 


SENSOR 


1,2,3,  and  4 IR 

5 , 6 , 7 , 8 , and  9 SLR 

10,11,  and  12  MTISLR 

13,] 4,  and  15  MTIFLR 

16,17,18,  and  19  ECM 

20  through  30  inclusive  PHOTO 

31,32,33,  and  34  FLIR 

35  Control  card  for  hard 

copy  output 

7.3.2  Outputs 

The  Executive  program  prepares  two  output  tapes: 

Output  Tape  1 describes,  on  a pass-target  basis,  the  detectability, 
and  localizability  statistics  for  each  of  the  sensors.  This 
tape  serves  as  input  to  the  Evaluation  Model  discussed  in  AIRS, 
Section  V.  The  second  tape  (Output  Tape  2)  contains  information 
computed  within  the  Executive  program  which  is  used  by  the  PI-0 
model  discussed  in  Section  V. 

In  addition,  if  the  user  so  desires,  a hard  copy  of 
the  information  recorded  on  Output  Tape  1 is  prepared. 

7.3.2. 1 Output  Tape  1 . This  tape,  also  called  the  Target/ 

Sensor  Output  Tape,  is  written  in  binary  form  and  is  organized 
on  a pass  basis.  Within  each  pass,  each  target  is  considered 
in  order  (up  to  a maximum  of  200  targets)  and  a group  of 
20  records  for  each  pass-target  combination  generated. 

The  20  records  for  each  combination  can  be  logically 
divided  into  three  groups: 


. Group  1 consists  of  the  first  15  records. 

Record  1 contains  general  header  information, 
true  target  data  (e.g.,  actual  time  of  beam 
passage),  and  various  system  parameters  (e.g., 
navigational  or  speed  errors).  This  information 
is  in  reality  the  statistics  provided  by  the 
Dummy  Sensor.  The  remaining  records  contain, 
for  each  of  the  four  levels  of  processing, 
statistics  on  the  14  other  sensors  - ordered  by 
sensor  number.  These  statistics  include  de- 
tectability, identif iability  , and  local i zabi 1 i ty , 
terrain  shadowing  probability  for  the  sensor 
and  target  being  considered,  etc. 

. Group  2 contains,  in  logical  records  16  through 
19,  the  combined  data  from  all  sensors  on  each 
level  of  processing  for  the  given  target.  This 
data  describes,  for  the  overall  system  of  sensors, 
how  effectively  the  target  was  detected,  identifie 
and  localized. 

. Group  3 consisting  of  record  20,  contains  the 
geodetic  probability  and  local i zabi 1 i ty  for  the 
target . 


Table  7-5  defines  these  records  in  greater  detail. 
Figure  7-1  shows  the  logical  grouping  of  records  on  the  tape. 

After  the  last  pass,  target,  and  sensor  have  been 
analyzed,  and  end-of-file  record,  which  does  not  contain  data, 
is  placed  at  the  end  of  the  tape. 
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7. 3. 2. 2 Output  Tape  2.  The  reader  is  referred  to  paragraphs 

3.6.2  of  AIRS,  Volume  I for  a complete  discussion  of  this 
tape . 

7. 3. 2. 3 Hard  Copy  Output . As  reported  in  paragraph  3.6.3  of 
AIRS,  Volume  I,  as  an  option  (if  computer  console  sense 
switch  HI  is  on)  a printout  of  Output  Tape  1 can  be  produced. 
The  output  is  easily  interpreted  since  all  variables  are 
labled.  This  output  does  not  print,  however,  words  35 
through  56  of  Croup  1 (Table  7-4)  as  data  contained  in 
these  words  will  be  processed  and  output  by  the  INTER 
program,  as  discussed  in  Section  VI  of  this  report. 
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latitude  and  longitude  are  to  be  punched  as  XXX.YXZZ,  -where  XXX  are  the  degrees,  IT  the  minutes 
nd  ZZ  the  seconds.  This  convention  applies  throughout. 

'ote:  If  target  type  I is  not  a radar,  card  may  be  left  blank. 

(!)  trust  be  specified  as  an  integer  between  -3  and  +3. 
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I JIR=1,  a diagnostic  printout  will  be  generate 
in  the  IR  subroutine.  If  JIR^l,  no  printout  is 
generated . 
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TABLE  7-3  INPUT  TO  EXECUTIVE  MODEL  (Cont. 
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TABLE  7-3  INPUT  TO  EXECUTIVE  MODEL  (Cont. 
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AL2(2)  lines  per  millimeter  for  sensor  1+6.  (Resolution 

taking  into  account  vibration,  noise,  turbulent  air  SF10. 
flow,  etc.) 
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TABLE  7-3  INPUT  TO  EXECUTIVE  MODEL  (Cont. 
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• • 

rH 

p 

a 

S-t 

c 

P 

o 

• • 

rH 

ex 

o 

f-> 

C 

ex 

o 

O 

rH 

rt 

V) 

• rt 

•H 

P 

o 

Pu 

M 

CO 

•rH 

M 

O 

Cl. 

M 

CO 

•H 

M 

4-> 

•H 

rH 

a5 

d 

rH 

M 

•rH 

to 

M 

MS 

CO 

CD 

M 

•H 

V) 

M 

Mi 

CO 

U 

■1J 

f-< 

O 

•M 

C 

•H 

C 

bO 

4-> 

rH 

3 

3 

o 

3 

O 

3 

+~> 

rH 

2 

3 

O 

3 

o 

3 

0 

rt 

CJ 

> 

Oj 

o 

o 

a 

o5 

r— 1 

r-l 

rt 

r-l 

H 

i — i 

05 

05 

rH 

1 — I 

rt 

r-l 

5h 

r-l 

<—* 

•rt 

> 

o 

Q 

u 

05 

u 

rg 

& 

£ 

P 

o 

X 

P 

P 

u 

£ 

ex 

O 

X> 

P 

ex 

o 

H 

<4-1 

o 

r-3 

to 

6 

<D  <D 
M > 

V ) O 
X-J 

to 

c 
o o 

> -H  -rf 
•r-t  M 
M p v5 
O H 
O O •> 

P CO  to 

<M  O •> 

W P CN1 


1/5  r-l  M 

rt  U C 
M O 
M C t/5 
COO 

o x p 
u s a 


d 

o x 

COO 
•rt  3 r— 1 
d O 
^ Ofl 
O M 
£ M C 
O C O 
P-rt  CO 

c 

o n o w 

C/5  M t/5  O 
•rt  M C 
O O X—l 
2 


d 

O X 
C O rt 

•rt  3 r-l 

d O 
h Otl 
O Si 
2:  M S-. 
O C O 

P*rt  CO 

c 

O CO  o 
Ui  M to 
• rt  M 
O O X 
2 


M d 

o o • •> 

CD  CO  CO 
M O rX 
CD  P O 

TJ  u o 

> H N 
rrt  rt  » 

rt  rH  r-l 

M X»— 1 
O M rt  to 

MM  r-l 

rrt  M O 
O M O > 
X3X  > O 
t-  O O Ml 


I CO 
M MS 
O O 
> O 

r-l  O iM 
rt 

C Xr-4 
O M 1-1 
•H  *rt  rt 

Mr-1  C4 

•rt  •»-(  O I 
d X>  O r-< 

c o > 

O M O CO 
U O r-l 

o d o 

O M O > 
X O C/3  O 
H d « -4 


to 

m 2; 
c 

CD  <W 
CO  1-H 

o 

>1 

CX  • to 

CD  O M 

>-  p c 
x o 

J5  M to 
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. -j.  - j.  -.4  \f-  V 


0)  ao 
U tO 

C -H  tO  to 
do  to  OO  P4 
0-r  bo  O 
do  >s  cO  O 

CH  4->  Oh  rH 
0-1  -H  GO 
dO  i— I rH 
•U  -H  CO  tH 
C .Q  CO 
•H  tO  • r 
■ JO  GO  S-i 
c o xj  o 
O Oh  O > 
<S  CXrH  o 


GO  r-l 

00  CO 

C +J  > 

CO  GO  ” 
Oh  to  4-> 

0-t  10 

+->  cO  cO  - 
c +->  >h  ^ 
CO  o 

•H  O +-»  O 
C/J  +-»  CO  r-l 

1 

+J 

U rH 
d>  rH 
4J  CO 
d)  CM 
t>  >1 
tO  O 

H c > Ifl 

CO  CO  O P<C 

era  o 

O tJ  o 

•H  X O rH 

+-»+->  to 

• H tH  CO  nd 
t)  H ^ 'H 
C -H  O rH 

o jo  J>  ri 

CJ  CO  CO  > 


C cO  tO  to 
d>  E o ^ 

Oh  00  O 

O >,(0  o 

4->  Oh  rH 


+-*  -H  CO  r— t 
to  JO  CO 
•r-t  CO  ' — 

i .n  <u  u 

C O pa  dO 

o Oh  o > 

21  (IH  O 


TO  CO 

r-l  C > tO 
CO  CO  O M 
C CO  O 
O tO  O 

•H  X O rH 
4~>  +->  00 

•h  -h  a rj 
to  i— I Oh  -h 

CH  OrJ 

O JD  > CO 
U CO  cO  > 


do 

to 

■H 

1 

Oh  l 

CO  CO 

o 

V) 

l/) 

Oh 

1 +->  1 Oh 

td 

o 

dO 

Oh  fd  Or  00 

o 

C > GO 

•H  o +-> 

>NtO  O 

P d)  c 
•h  to^-- 
i — i c0 
•H  Oh  CO 
X>  d>  ±L 

to  > O 

POO 


+->  x>  rH 
d)  OH 
Q H O 


• H -H  tO 
3 «H 
r CT  tO 
MO  fi 

co  o 

•H  Oh  t_> 

o 


rH  00  O 
JO  U-3 

c 

d o p c 

•H  Oh  H G 
CO  O O *r- 

i-  i+h  e +j 


CM  -tO 

o to 

t3  C0-H 
(0  rtH 
cO  Ih  B 

ea  o > 
• > 

+->  CO  rH 
to  rH 
tO  to  CO 
Oh  tO 
4->  3 Oh 

d o o 

O rH  > 


Oh 

O I 

to  d) 

C Oh 
d)  d 
to  d> 

Oh 

10 

+->  X 

e 

dO  4-1 
V)  I— i 
dO 
Oh 

Dr  • tO 

do  do  4-> 
Oh  EX  d 
>.  CO 
2d  4->  CO 
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TABLE  6-3  DATA  ARRAY  FOR  EACH  PASS-TARGET  COMBINATION  (Cont. 
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APPENDIX  A. 

^ PROOF  OF  ALGORITHM  FOR  DI-TERMINI  NG 
PROBABILITY  OF  REALIZING  k OUT  OF  N'h  EVENTS 


V 

V 
>< 
N 


r. 


We  first  find  the  probability  of  exactly  k targets 
identified  as  type  T^.  Let  be  the  number  of 

t%  of  type  T,  in  the  complex.  Let  E.  = the  event 

. th  n J 

i\yc  j target  of  type  T^  is  identified  where 

- *V 


V 


Let 
occur  (i.c 


E.  . = the  event  that  exactly 
K » J 


■ em 


k events 
= the  event  that  exactly  k 


J <?  out  of  the  first  j targets  of  type  are  iden*We«A 

► . Thus  we  are  seeking  first  P(Ej*t^  ),  the  probab^*^ 


of  exactly  k 


targets  of  type  T^ 


being  identified 


all  N targets  of  type  T^.  We  find 
by  iteration  (induction)  over  j. 


P(Ek  •)  for 

K » J 


A-  1 


t w * . • « * V f 


4 


or 


V V. 
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In  other  words,  we  first  find  P (Ej,  q)  for 
k = 0>N}j-  Then  we  calculate  P(Ej,  for  k = 0»N^  fron 

P(Ej,  q)  for  k = 0,N^  using  the  general  technique  described 
below.  Continuing  in  this  manner,  we  finally  calculate 

p(Ek,Nh)  for  k = °»Nh  £rom  pCEk,Nh-l)  for  k = °»Nh‘ 


Clearly  Eq  q always  happens,  since  exactly  0 
targets  of  type  will  be  identified  out  of  the  first  0 

targets  of  type  T^.  Similarly,  E^  q k f 0 never  happens 
since  exactly  k targets  out  of  the  first  0 targets  of 
type  T,  cannot  be  identified  if  k * 0 . Therefore, 


cannot  be  identified  if  k *■  0 . 


pCEk,cP  = 


1 if  k = 0 


0 if  k * 0 


(A.  1) 


So  we  know  P(E^  q)  for  k = 0,N^.  Now  suppose  we  know 
P(E.  .)  for  k = 0,N,  . We  want  to  calculate  P(E,  . 
for  k = OjN^.  First  we  find  PfE^  j+i^  £or  k = 0. 

E0  j+1  happens  only  when  Eq  ■ happens  and  Ej+i  does  not 
happen.  So  p(E0j  + 1)  = P^EkJ^1  " P^Ej+l^’  Next  we 
calculate  P(Ej.  j+i^  for  k > 0*  We  distinguish  between 
two  mutually  exclusive  and  exhaustive  cases  of  Ej.  ^ + ^ 


(1)  Ej+x  occurs  and  occurs 

(2)  E . , does  not  occur  and  E.  • occurs 

J +i-  K > J 

The  probability  of  case  (1)  is  P (E . . ) P(E.  . since 

j + j.  k-  i , 3 

Ej+^  and  E^.^  • are  independent.  Likewise,  the  probability 
of  case  (2)  is  1-P(E.  ,)  P(E,  -).  In  addition,  P(E,  - .) 

j + -*■  K y J ^ > J -» 

is  the  probability  of  case  (1)  plus  probability  of  case  (2). 
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Therefore,  we  conclude: 


P<Bk,j,l)  - p(Ej  + iJ  p(Ek-l,j)  * , for  k> 


Thus  we  now  have  the  general  technique  to  find  P (E^  . ) 

for  k = 0,Nh  in  terms  of  P(Ek  ^ ) for  k = 0,Nh-  Now’3 

by  Nh  successive  applications  of  this  technique  we  find 

w ) f°r  k = 0,N\  . This  is  our  desired  result. 
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APPENDIX  B 

«£,  PROBABILITY  THAT  AT  LEAST  ONE  OF  SEVERAL  DECISION 
FUNCTIONS  IS  SATISFIED 


This  appendix  describes  in  detail  the  algorithm  used  ytt 
determine  the  probability  of  observing  the  target  complex. 

the  decision  rule  is  outlined.  Then  the  algorithm  is  • 

l-QIDAiiL 


The  decision  rule  has  associated  with  it  target  types 

T?  > • ■ ,T  and  decision  functions  D.  , D~,...,D  . The  rule  is 
c g i 4 , n 

isfied  (i.e.}  the  complex  is  observed)  if  at  least  one  of 
decision  functions  is  satisfied.  Each  decision  function 
sists  of  a set  of  g numbers,  each  corresponding  to  a target 
e. - Each  number  specifies  the  minimum  number,  i,  of  targets 
,the  corresponding  type  which  must  be  identified  for  the 
lision  function  to  be  satisfied. 


J ( ) Jv’e  are  given  the  probability  of  identifying  at  least 
|=0, 1,2...)  targets  of  each  type  T1,T2<..T  , and  wish  to 
jrmine  the  probability  of  observing  the  complex.  Obviously, 
limit  the  decision  rule  to  the  first  target  type,  T^ , the 
(ability  of  observing  the  complex  would  be  trivial  to 
vite.  In  that  case,  letting  m be  the  least  minimum  (for  the 
; type,  Tj)  of  all  the  decision  functions,  we  have  immediately 
■ the  probability  of  identifying  at  least  m is  the  probability 
^serving  the  complex.  We  could  limit  the  decision  rule 
jily  to  the  first  type  T, , but  still  further  to  some  subset 

f *■ 

V set  of  decision  functions.  In  this  case  the  probability 
'icrving  the  complex  is  obtained  just  as  trivially.  Since 
»*se  with  the  rule  limited  to  T^  is  so  simple  it  is  natural 
consider  the  case  with  two  target  types  T^  and  T2- 
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It  is  shown  later  in  this  appendix  that  the  probability  of 
observing  the  complex  with  the  rule  limited  to  and  1’2  can 
be  obtained  if  one  knows  the  probability  of  observing  the 
complex  with  the  rule  limited  to  a subset  of  the  decision 
functions  which  in  turn  are  limited  to  . In  general,  it 
is  shown  that,  if  we  can  calculate  the  probability  of 
observing  the  complex  with  the  decision  rule  limited  to  a 
subset  of  the  decision  functions  which  in  turn  are  limited 
to  types  T^,T2,...,T  , then  we  can  calculate  the  probability 
of  observing  the  complex  with  the  decision  rule  limited  to 
some  subset  of  the  decision  functions  which  are  in  turn 
limited  to  types  ,T2 , . . . ,T  ^ • Continuing  with  this  pro- 
cedure we  obtain  the  probability  of  observing  the  complex 
with  the  decision  rule  limited  to  some  subset  of  the  decision 

i 

functions  which  are  in  turn  limited  to  types  T^,T2,...Tg. 

If  we  let  the  subset  of  decision  functions  include  all  the 
decision  functions,  we  have  not  limited  the  rule  at  all. 

Hence  we  would  then  have  calculated  the  true  probability  of 
observing  the  complex. 
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To  explain  the  above  procedure  we  define  some  useful 
terms.  Let  Ev  . = the  event  of  satisfying  decision  function 

D,  ' s minimum  for  target  type  T. . Notice,  E,  . is  merely  the 
event  that  at  least  m (the  minimum  for  target  type  in  Dj.) 
targets  of  type  T\  are  identified.  Furthermore,  for  all 
m=0,l,2...  we  know  the  probability  of  identifying  at  least 

m targets  of  type  T\.  # Hence  we  know  P (Ej,  , the  probability 

that  E,  • occurs.  Also  E,_  . is  the  event  that 

k 1 1 i<i  K > J 

decision  function  is  satisfied  for  all  target  types 

T, ,T-, . . . ,T. . In  particular  .Q  E,  . is  the  event  that 
1 ’ 2 l F J - g k , j ,n  p | 

decision  function  Dj.  is  satisfied,.  Consequently  j<  o 

E,  . is  the  event  of  satisfying  at  least  one  of  the  dccison 
k > I 

functions,  i.e.  the  event  of  observing  the  complex.  Similarly, 


0 

k=l 


n 


E,  . is  the  event  of  satisfying  the  decision 
^ > J 

rule  when  it  is  limited  to  target  types  ,T? , . . . , . 


j < i 


B-2 


A 


m 


Now,  let  SC  |l,2,...,n|  so  S represents  a subset 
(the  set  of  all  decision  functions  with  index  in  S)  of  |D,,D2, 
...,Dn|  the  set  of  decision  functions.  Then  '<i 


the  set  of  decision  functions.  Then 


Eu  ,.  is  the  event  of  satisfying  the  rule  when  it  is  limited  to 
K > J 

the  subset  of  decision  functions  represented  by  S,  which  are 

in  turn  limited  to  types  T, ,...,T..  This  event  is  referred 
lx — 1 i I | n 

to  frequently,  so  we  devise  a new  symbol:  FQ  . = , U 

E,  ..  By  definition,  F c.  0 » is  the  event  of  observing 

k,j  7 \ 1,2, . . . ,n } ,g  6 

the  complex.  So  we  want  to  calculate  P(I;  i > ) 

|t,...,nj  ,g 

This  is  done  by  iteration  (induction)  over  the  second  index. 

That  is,  each  step  of  the  iteration  increments  the  second 
index*,  thus,  as  mentioned  earlier,  incrementing  the  number  of 
types  to  which  we  restrict  the  rule. 
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To  perform  the  iteration  we  must  first  calculate 
P(Fo  • ) w'here  S is  a non-empty  subset  of  {l,...,n{  . (0  is 

the  empty  set  so  we  write:  0 0 S C {l,...,n|  ). 


For  each  k«  s we  have  a minimum  number  of  targets 
required  to  be  identified  as  type  T^  to  satisfy  Dj, . he  may 
order  the  decision  functions  in  S by  the  minimum  number  of 
targets  they  require  to  be  identified  as  type  Tj . Let 
be  the  decision  function  requiring  the  least  number  (that  is 
the  most  easily  satisfied)  of  type  T^ ; be  the  one 

requiring  the  second  least  number  of  type  T^ ; etc.  We  notice 
that  if  Fg  j occurs  Dqj,^  is  satisfied  for  type  , because 
if  any  decision  function  with  indice  in  S is  satisfied  for 
type  Tj , so  is  b0RUJ  . Also,  if  is  satisfied  for  type  T^ 

Fg  ^ occurs  by  definition.  Therefore,  P(Fg  ^)  = P(Pqrd1  i ^ 
for  all  non-empty  S C |l,2,...,n|  . However,  as  explained 

earlier  we  know  P(bQRD  ) and  therefore  we  know  P(Fg  1). 

Now  wo  describe  the  general  technique  used  in  the 

iteration  procedure.  Suppose  we  know  P(F<.  ) for  all  non-empty 

o , r 
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sc  { 1 , 2 nj  and  want  to  find  P(FS  where  0 f SC|l,2, 

We  first  reorder  the  decision  functions  for  Tr+1  as  we  did 
above  for  T^.  Let  the  number  of  elements  in  S be  2.  Then 
for  calculation  purposes  we  break  Fs  r+1  into  the  following 
mutually  exclusive  and  exhaustive  cases: 


1.  F 


| 0RD1  } , r 


"OilDl , r + 1 


"0RD2 , r+1 


{ 0RD1 , 0RD2  J , r 


J0RD2  ,r+l 


J0RD3 , r+1 


occurs  and 


occurs  and 


does  not  occur 


occurs  and 


occurs  and 


does  not  occur 
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q-1.  F 


{ 0RD1 , 0RD2  , . . . ,ORD(q-l)|  ,r 


occurs  and 


JORD(q-l) ,r  + l 


JORDq , r+1 


occurs • and 


does  not  occur 


^ ’ P | 0RD1 , 0RD2 , . . . ,ORDq  j , r occurs  and 


"ORDq , r+1 


occurs 


| « 

V ■>  • 

V;^ 

vj.  *>% 

Sr'Vv 

:-rvVv 


1.  The  probability  of  case  1 


P(F  { 0RD1  | , r^  PfIZORDl,r  + l)  ' P (P0RD2 , r+1 ^ 


VWA'-A 


^ ' %.  -.  T 


2.  The  probability  of  case  2 = 

P(F  | 0RD1,0RI)2  | , r ^ [P(E0RD2, 


r + 1^  ' P(E0RD3 


q-1.  The  probability  of  case  q-1  = 


» r+1 '1 


P(F  |0RD1,. . .ORD(q-l)  } , r}  [P (EORD(q-l) ,r+l3 


P(EORDq,  r+1 
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q.  The  probability  of  case  q = 


E 


P(F  { 0RD1 , 0RD2 , . . . ORDq  | , r^ P (EORDq , r+1^ 

All  of  the  values  on  the  right  of  equations  1 through 

q are  known  so  we  know  the  probability  of  each  case.  Also 

P(FS  r+1)  = probability  of  case  1 +...+  probability  of  case  q, 

since  the  cases  are  mutually  exclusive  and  exhaustive.  Thus, 

we  now  have  the  general  technique  for  finding  P(FS  +1) 

tvhere  0 f S C{l,...,n{  given  PCF^  ) for  non-empty  S cjl,...,i 

By  g - 1 successive  applications  of  this  technique  to  P(FQ  ) 

for  non-empty  S Cjl, . . . ,n  | we  obtain  P(Fq  ) for  non-empty 

S Cjl,...,n|  . In  particular,  we  obtain  P(F  <.  . ), 

j 1 , . . . ,n  { ,g 

the  probability  of  observing  the  complex. 
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To  be  certain  of  the  validity  of  our  result  we  must 
prove  our  assertion  that  the  q cases  arc  indeed  mutually 
exclusive  and  exhaustive.  To  prove  the  q cases  mutually 
exclusive,  consider  two  different  cases  J and  K.  Without  loss 
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of  generality  we  may  assume  J < K.  Case  J requires  that 
EORD(J+]),r+l  does  not  occur>  which  implies  that  ^0RDKjr+1 
do^.'S  not  occur.  (This  is  clear  from  the  following  argument. 

E01!D(J*I),r*l  doos  not  °ccur  moans  that  D0RD(J*1)  is  not 
satisfied  for  target  type  Tr+1,  but  D0RDfJ+1)  is  satisfied  for 
target  type  Tr+1,  whenever  D0RM  is  because  J+l  < K.  So 
EORDK  raust  not:  i,e  satisfied  for  target  type  Tr+^;  that  is 
EORDK,r+l  does  not  occur.)  However,  case  K requires  E0RDR  + , 
to  occur.  Therefore,  case  J and  case  K cannot  both  happen  at 
once.  Consequently,  the  q cases  are  mutually  exclusive. 


To  show  the  q cases  are  exhaustive,  suppose  that 
Fs  r+2  has  occurred.  Then  there  is  some  K such  that  ORDK * S 

and  DORDK  is  satisfied  for  Ti » • • • >Tr+i • Let  DQRDJ  be  the 
last  decision  function  in  (D0R1)1>  D0RD2>...,  D0RDq)  that  is 
satisfied  for  target  type  T j.  Then  case  J is  satisfied  as 
seen  by  the  following  inspection  technique.  F (nnni  onn.  \ 

j UKU  I j • a • y Ul\Dj  J ) 2 

must  occur  since  K < J (because  DqR^r  is  satisfied)  and 
DqruK  satisfied  for  T^,...,T  ^ by  assumption.  Also 

EORDJ,r+l  occurs  since  DqRRj  is  satisfied  for  target  type 


*r+l  * Similarly,  e0rd(J+1) ,r+l  UUCb  UUL  uccur  5iace  uORDJ 
was  chosen  as  the  last  decision  function  to  be  satisfied  by 

target  type  Tr+1 ' Thus,  the  q cases  are  also  exhaustive. 


does  not  occur  since  D, 


Therefore,  our  method  is  indeed  valid  and  by  q - 1 
successive  applications  we  obtain  the  probability  of  observing 
the  complex. 


